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PREFACE 


This is an interim report covering the research completed under con- 
tract grant number NAS8- 29-100 titled "Water Resources Planning for Rivers 
Draining Into Mobile Bay. Part II: Non- Conservative Species Transport 

Models." This report covers the period January 1 to December 31, 1974, and 
serves as Mr. Hua-An Liu’s M.S. thesis which will be used as partial fulfill 
ment of the requirements for that degree at The University of Alabama. 

A third report, "Part III*. Application of Developed Models to User 
Needs" will be issued at the termination of the next grant period. 
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ABSTRACT 


The purpose of "this research effort is to expand the 
mathematical modeling capabilities of the hydrodynamic and salinity 
models of Hill and April to include a description of non- conservative 
species transport in the Mobile Bay system. In so doing, the 
knowledge gained provides a clear insight into the effect that rivers 
draining into the bay have on water quality conditions. 

Total coliform group bacteria were selected because of their 
relationship to commercial fishing ventures within bay waters. This 
item was also chosen on the basis of data availability sufficient 
for model calibration and verification. Results are presented as 
monthly average distributions corresponding to the data base used. 

In addition to the above, a parametric study was also conducted. 
In this study river flow rates, wind conditions and bay system tem- 
peratures were investigated to determine their influence on the total 
coliform concentration patterns. Of these factors temperature and 
river flow rate had a pronounced effect on the concentration profiles, 
while wind conditions showed only slight effects. Shifts in concen- 
tration profiles as much as 8 kilometers were observed in extreme 
cases. 

The effect of changing total coliform group loading concen- 
trations at constant river flow rates and temperature was also inves- 
tigated. As expected these loading changes had an appreciable 
influence on total coliform distribution within Mobile Bay, 

iv 



Utilization of the Non-conservative Species Transport Model 
to predict trend "behavior in the Mobile Bay system is demonstrated* 
Continuing efforts to improve the data collection programs in support 
of mathematical modeling are encouraged to increase the utility and 
predictive capabilities of the models. 
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CHAPTER I 


INTRODUCTION 

Sewage, industrial waste disposals, and storm water overflows 
discharged into Mobile River and surrounding creeks from the Mobile 
metropolitan area, and excessive concentrations of bacteria in the 
Mobile River, result in the pollution of Mobile Bay. A location map 
showing these sources of waste is shown in Figure 1-1. One method 
for expressing the bacterial content of these waters is to determine 
the total coliform bacteria group count which gives an indication of 
the disease carrying bacteria or pathogenic content in the water. 
Because of this pollution, Alabama, under state laws and the regu- 
lations of the State Board of Health, closes the bay to oyster 
harvesting as a safeguard to human health. The criterion on which 
closing the bay is based is either a total coliform concentration in 
excess of 70 parts per 100 ml at specific locations adjacent to oyster 
reefs, or whenever the concentrations of 10% of all samples collected 
are in excess of 230 parts per 100 ml. ^*6) These samples are 
obtained in the field and analyzed in the State Laboratories at 
Montgomery, Alabama. In current years, this policy has led to the 
permanent closing of the upper third of the bay, the intermittent 
closing of the middle third, and closing of the lower third during 
extremely high pollution periods. These sections of the bay are more 
clearly defined in Figure 1-2. From 195^ to 196? bay closures 
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resulted in annual losses of approximately a quarter of a million 
dollars from oyster harvesting alone. 

It becomes apparent, from the economic considerations 
associated with the maintenance of safe oyster harvesting conditions, 
that a rapid predictive method, supplemented with spot analytical 
support, could result in substantial savings of time and effort. 
Furthermore, the method could provide a way of determining the effects 
system variables, such as river flow rates, runoff, degree of waste 
treatment, and expansion of waste treatment facilities, have on the 
coliform distribution in the bay. This technique could also provide 
clues as to ways in which these most serious upsets to the bay could 
be abated. 

This study provides such a method which has as a basis the 
application of conservation of mass and species equations subject to 
the bay ecosystem constraints. For this purpose, a two dimensional 
(surface), non- conservative species transport model is developed for 
Mobile Bay. The model is solved with a finite difference method and 
implemented by computer solution using a UNIVAG 1110 system. The 
hydrodynamic model for Mobile Bay developed by Hill and April is 
used to provide basic current and dispersion coefficient data required 
by the non- conservative species transport model. The resultant 
package, referred to as the Non- conservative Species Transport Model 
(NCSTM) is verified with available total coliform bacteria data 
obtained from the State Department of Health. Extension of the 
NCSTM can be made to analyze the biological oxygen demand (BOD) and 



5 


dissolved oxygen (DO) levels within the "bay. Model verification for 
BOD and DO are deferred until field data become available. 

Parametric studies are included to determine the effect that 
system variables such as wind speed and direction, river flow rates 
and temperature have on the coliform concentration distributions 
within the bay. Based on these studies conclusions are drawn which 
indicate the conditions most conducive to pollution flushing and 
dispersion in the Mobile Bay ecosystem. 



CHAPTER II 


BACKGROUND 

Due to the complex nature of estuarine systems it very often 
is not feasible or practical to study the behavior of the systems 
by field data analysis. Many sampling stations must be monitored In 
a way to determine meaningful results about what one part of the 
system is doing relative to another. These so called synoptic sam- 
pling plans require a great number of research vessels and man hours 
to obtain accurate and precise data to determine the real behavior of 
the system. Additionally, during periods in which bad weather occurs, 
the system data collection plan is often inoperable. Mathematical 
and physical modeling of these systems have been demonstrated to be 
reasonable methods to circumvent these problems. 

2.1 Model Concept 

A model is, in short, a representation of the real system. 
Various models have been used to study the hydrodynamic behavior and 
water quality conditions of streams and estuaries. An acceptable 
model is one in which specific responses caused by variations in 
system parameters can be reasonably and accurately described. In 
order to show acceptability there are two phases which must be de- 
monstrated when models are utilized. These are the calibration and 
the verification phases. 
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Some characteristics of the real system may not be sufficiently 
understood and some empirical equations may be required to correlate 
the resulting behavior. These correlations would substantially 
depend on the specific real system and may vary from system to system. 
Before the model can be verified, it is necessary to find the set of 
correlations which best describes these characteristics for the 
specific system under consideration. This is called model calibration. 
After the model is calibrated, the correlations are fixed to perform 
the verification of the model. The use of a model to successfully 
predict what would happen in the real system due to variation in 
system parameters for a given period results in verification. This 
phase requires the availability of sound data to show that model 
predicted results are in fact duplicative of system behavior. Failure 
to achieve comparative results during this phase of the study could 
result in either recalibration of the model or collection of field 
data more representative of the real system behavior, or the develop- 
ment of a new model. Statistical analysis during this phase of the 
study is essential. 

Once verified, sensitivity of the model predicted results can 
be studied by a parametric investigation. In this phase of the 
project, system variables thought to be important can be varied 
individually with measurement of the response in the objective 
function. Significant changes in the objective function ( in this 
study the concentration of non- conservative species ) for each 
perturbation of the independent variable are then a measure of its 


sensitivity. 
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These phases of the study are intended to establish confidence 
in the model predicted results. 

2.2 Modeling Estuarine Systems 

Generically, models used to describe estuarine systems can be 
divided into two types, physical models and mathematical models. 

A physical model is a scaled imitation of the real system. 

There is a physical model for Mobile Bay at the Water Experimental 
Station of the Corps of Engineers at Vicksburg, Mississippi cons- 
tructed in 19?3 at a cost of approximately $1,000,000. It has been 
successfully used to reproduce tidal and current conditions and to 
simulate dispersion effects with dye tracer release experiments which 
in turn provide useful information about mass transfer rates in the 
bay. Characteristics of the physical modeling have been discussed 
by Masch(^®)« 

A mathematical model is a functional representation of the 
real system, i.e. a set of partial differential equations describing 
the system under study and the associated assumptions and constraints 
that apply to its formulation. Mathematical models can be divided into 
analog models and digital models according to the type of computing 
facilities used to implement the numerical solution to the partial 
differential equations. With the development of high speed digital 
computers, mathematical models using finite difference methods to 
solve the partial differential equations have become widely accepted. 
The model used for the study of Mobile Bay is a mathematical model 
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implemented, by a high speed digital computer. 

Many mathematical models for estuarine water quality have 
been developed for various systems (25) , These models are further 
classified in terms of the spatial and temporal conditions over which 
they are designed to perform. These include one-, two- and three- 
dimensional steady and transient models. The application of the 
specific model to be used is dependent on the system geometry, hy- 
drology, and the time frame for which information is desired. 

One-dimensional Models 

The transient species continuity equation for one dimensional 
systems can be written as 

+ (2-1) 

where c = concentration of the water quality species along the 
direction of stream flow 
t = time 

u = velocity of stream over the cross-section of flow 

E = dispersion coefficient 

x = distance in the direction of flow 

Sj_ = sources or sinks of the water quality species 

For narrow waterways where cross-sectional variations in 
physical and water quality parameters are negligible, such as creeks, 
rivers, and narrow estuaries, the one-dimensional model is justified. 
Again, due to the complexities in the physical systems, complete 
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analytical solutions are not always possible. Two approaches, i.e. 
the continuous solution approach and the finite section approach, 
have been utilized in solving one-dimensional problems in estuaries. 

In the continuous solution approach, it is necessary to divide 
the system into a number of individual sections or subsystems, each 
of which is characterized by physical and hydraulic parameters. 

Sections are joined by related concentration and flux terms. Analy- 
tical solution of the one -dimensional equation (2-1 ) may then be 
obtained for each section \ they are then summed up to give the overall 
solution for the system. This approach was adopted in the East River 
Model ^ . 

With the finite difference approach the differential equation 
is replaced with a difference equation and the system is divided into 
a number of sections with the assumption of complete mixing in each 
section. Matrix inversion or relaxation is then used to obtain 
solutions. This approach was used by the Thames Estuary Model, the 
Delaware Estuary Model, the Potomac Estuary Model, and the San Eran- 
cisco Bay System Model which consists of San Francisco Bay, San Pablo 
Bay, Suisun Bay, and Sacramento -San Joaquin Delta. 

Two- and Three-dimensional Models 

The two-dimensional transient species model equation can be 
written as 

= — (u ~~ + v -fy-) + -"(Z-z) 

where x and y are the two dimensions over which c is varying. 



For water bodies which are well mixed in one of the three 
dimensions, the use of two-dimensional models which describe the 
variation of conditions in the second and third dimensions, is 
justified. Most wide, shallow bay systems which are vertically well 
mixed fall into this category. Several examples of modeling studies 
involving bays in the United States are summarized below. 


(i) San Francisco Bay System Model - 


This model was principally 


developed by Water Resources Engineers, Inc. The basic approach 
was to represent the estuary with a network of uniform channels 
interconnected at junctions or nodes . This allowed a one-dimen- 
sional treatment of a two dimensional system. It has been effec- 
tively and extensively applied to the San Francisco Bay system. 
The computational experience thus developed was then utilized 
in modeling Sydney Harbor, San Diego Bay and the Columbia River, 
etc. Verification was made on salinity data & 5 ) , 


(2) Hillsborough Bay Model The finite difference approach as 

applied to the one -dimensional model was extended to two dimen- 
sions in Hillsborough Bay which is a natural arm of Tampa Bay, 
Florida. The bay was horizontally segmented such that in each 
segment the depth is approximately uniform. Verification was 
made on bay salinity and investigations were performed to study 
the effect of diverting the Hillsborough River as a means of 
smoothing the fluctuations of salinity in the bay. 


( 3 ) Galveston Bay Model This is a time-dependent two dimen- 


sional model using finite difference methods to solve the model 
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equations. The "bay system as well as the Houston Ship Channel 
were segmentized into squares of uniform size. The model was 
extensively verified using salinity, BOD, and DO collected from 
-water quality stations maintained in the system. 

( 4 ) Louisiana Coastal Marsh Area Model Two-dimensional, time- 

dependent models were developed which predict the velocity 
profiles, tidal fluctuations, and temperature and salinity pro- 
files for the Barataria Bay region of coastal Louisiana. An 
alternating-direction implicit, finite difference method was used 
to solve the differential equations numerically. Results were 
reported for the dynamics of tidal fluctuations, velocity profiles, 
and salinity and temperature distributions for conditions encoun- 
tered in May 1 970 ^ ^ ) . 

(5) Mobile Bay Model This study covered the hydrodynamics and 

salinity of Mobile Bay, Alabama. It accurately predicted time- 
varying tidal heights, current patterns, and averaged salt con- 
centration distributions of the Bay. A two-dimensional finite 
difference method was used to approach the explicit solution to 
the model equation. A salinity wedge was used in the lower reach 
of the bay to simulate the saline water intrusion without going 
to a much more costly three-dimensional computational scheme. 

In some cases, in addition to changes in the horizontal dimen- 
sions, the flow may be highly stratified in the vertical direction, 
causing significant gradients in the depth direction. This pheno- 
menon is frequently observed at locations near waste outfalls. A 
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three-dimensional model would find particular application under such 
conditions . 

Any water quality model relies on its hydrodynamic counterpart 
for hydrodynamic parameters such as current velocities and dispersion 
coefficients. The basis for the Non -conservative Species Transport 
Model presented in this study is the Hydrodynamic model of Mobile Bay 
developed by Hill and April 

2.3 Non-conservative Species vs. Conservative Species 

The term "non- conservative species" is used to refer to the 
materials dissolved in the estuarine water in which the concentrations 
are subjected to rather rapid and appreciable changes. These changes 
are caused by various mechanisms of generation (source) and dissipa- 
tion (sink) depending on the characteristics of the species itself, 
the physical environments to which they are exposed, and other aquatic 
ecosystems with which they are interacting. For example, most water 
quality entities of great concern to us, such as total coliform, BOD, 
and DO, are very sensitive to physical, biological or chemical upsets 
which result in changes in the environmental balance within the sys- 
tem, They are essentially non- conservative in nature, and are gene- 
rally named "non- conservative species" in water quality studies. 

The "conservative species" label, on the other hand, is used 
to refer to materials dissolved in estuarine water in which the con- 
centrations are rather stable as compared to non- conservative species. 
Most of these species are not chemically or biologically reactive 
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substances. For example, salinity concentration is affected by 
freshwater discharge, rainfall, evaporation, and sea water intrusion, 
instead of any appreciable biochemical or chemical effects. However, 
in some cases salinity may also be considered a non- conservative 
species. A salinity model (conservative) for Mobile Bay has been 
developed by Hill and April 

2.4 Non- conservative Species Modeling in Bay System Analysis 

It has been estimated that approximately one third of the 
total population of the United States, or 40 of the 110 Standard 
Metropolitan Statistical Areas are located on estuaries (l?) . The 
vulnerability of estuarine systems to human influence has been 
demonstrated in recent years by observed upsets. Methods for the 
abatement of pollution of these delicate systems are being sought 
with increased intensity. With the advance of technology and the 
rapid growth in population, people are making much more use of the 
natural environment and at the same time dispose much more waste into 
it. Estuarine systems, which have long been depended on for their 
ability to assimilate a variety of wastes, are now becoming the first 
victims. Unlike the Olympic National Forest in the State of Washing- 
ton * w hich is known for its ecological stability, estuarine sys- 
tems are unstable, and subject to an increasing number of man-made 
and natural disturbances. 

While digesting the waste input from various sources, estuarine 
systems have to maintain their own natural balances. When changes 



are occurring gradually, estuarine systems can adjust to them quite 
well. However, present day upsets are occurring over extremely short 
periods, which overburden the evolutional process or homeostatic 
ability of the systems. These processes do not always have time to 
optimally operate, and the stability of the systems becomes critical. 
This threatens the existence of the ecosystems within the estuaries 
and seriously reduces the ability of the estuaries to provide people 
with those resources taken for granted for such a long time. This 
in turn affects the quality of life of the entire population. The 
yearly closing of Mobile Bay to oyster harvesting or the elimination 
of recreational activities in Lake Pontchartrain in Louisiana or the 
permanent restriction of waterways to navigation status only are just 
a few examples of loss of natural resources utilisation. Most upsets 
in these waters result from excessive waste disposal from municipal, 
industrial or agrarian sources. 

Water serves as a good medium for disease- carrying organisms. 
The bacteria of typhoid fever, cholera, and dysentery are all water 
borne pathogenes. It is assumed that the number of disease- carrying 
microorganisms in water is proportional to the total number of micro- 
organisms. Due to the variety of microorganisms, it is impossible 
to perform quantitative tests determining all the species. The total 
coliform bacteria group count, which is a count of the total bacteria 
content, therefore, becomes an indication of the disease -carrying 
bacteria, or the pathogenes within the water system. A high patho- 
gene content renders water hazardous to the persons using the estuary 
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for havesting, recreation, and even navigation. The total coliform 
concentration standard for shellfish harvesting in coastal and marine 
water adopted by the State of Alabama is "not to exceed a median MPN 
(most probable number) of 70/100 ml and not more than 10% of the 
samples shall ordinarily exceed an MPN of 230/100 ml for a 5 - tube 
decimal dilution," which is consistent with standards used by the 
National Shellfish Sanitation Program as well as some other states^) 
When these criteria are exceeded, the bay waters are declared hazar- 
dous to health and are closed to public use. This form of pollution 
is often seasonal, occurring during periods of heavy rainfall and 
correspondingly high runoff rates. 

Another kind of pollution results when discharges of organic 
materials occur. These organic materials can serve as nutrients for 
microorganisms. These organisms digest the wastes with the excretion 
of more elementary type materials which can serve as food to be ab- 
sorbed by phytoplankton and plants within the system. In these di- 
gesting processes oxygen is consumed. Therefore, when a sudden 
excessive amount of nutrients is introduced, the oxygen content may 
rapidly decrease to a very low level or even entirely vanish, because 
the reaeration mechanisms are not able to keep pace with the oxygen 
consumption rates. This total depletion of oxygen, although lasting 
only a short period of time, often results in fish kills. The "Jubi- 
lee" recorded in the northeastern coast of Mobile Bay v ' ' and some 
other parts of the Gulf Coast areas, are examples of this phenomenon. 
Under such circumstances, the organic materials introduced as wastes 
are no longer nutrients, but are instead pollutants. The control of 
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such -waste materials is predicated on the sound knowledge of the 
system behavior including those hydrodynamic , biological and chemical 
processes which describe its assimilating capacity. To analyse such 
behavior, description of those species which make up the system are 
essential. This investigation is directed at the development of a 
rapid, accurate, predictive method for describing non-conservative 
species transport patterns in Mobile Bay. 



CHAPTER III 


DERIVATION OF THE NON -CONSERVATIVE 
SPECIES TRANSPORT MODEL FOR MOBILE BAY 

The differential equations used in the Non -conservative Species 
Transport Model (abbreviated as NCSTM) are derived in this chapter. 

The general differential equation developed is modified according to 
spatial and temporal simplifications, and through characteristic 
constraints of the real system. Numerical form of the model equation 
is then presented together with the solution procedures. 

3*1 The Physical Setting 

Mobile Bay is approximately 49 km. (31 miles) long and has an 
area of 1070 km.^ (419 square miles)^^^. It has a ship channel which 
has a total length of 3^.5 miles and is 40 feet deep and 400 feet wide. 
The channel runs through the left half of the bay from the Main Pass 
at the Gulf of Mexico in the south to the Mobile River in the North. 

An intercoastal waterway, which is 12 feet deep and 200 feet wide, runs 
from west to east from Grant's Pass between Little Dauphin Island and 
Cedar Point toward the lower right comer (Bon Secour) of the Bay. 
Except for the ship channel and the Intercoastal Waterway, the Bay is 
shallow with a flat bottom. The average depth is 9*81 feet at mean 
low tide. Six rivers drain into Mobile Bay from its perimeter (see 
Fig. 1.2). Naming them in a counterclockwise manner beginning in the 
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northwest these rivers are the Motile, Dog, Fowl, Bon Secour, ,Fisn and 
Tensaw Rivers. The Mobile and Tensaw Rivers are the largest of the 
six, with average combined volumetric rate of discharge of 59>0OO cfs. 
The Mobile River perennially discharges large amounts of highly conta- 
minated waters and is considered the main source of pollution loading 
to the bay. Dog River, located near the Mobile River to the southwest, 
may also contribute substantially to the pollutant concentration in the 
bay. 


Average atmospheric temperatures over the Mobile Bay area have 
been accumulated by the Weather Bureau of the U. S. Department of 
Commerce (^7 ) . Monthly averages range from approximately $0 °F in the 
cold months to the low eighties in the warm months. Wind speeds and 
directions over the bay are also included in the climatological data 
collected by the Weather Bureau(^), Monthly averages range from 
approximately f mph. to 13 mph. for the period January to August, 1962. 


3*2 Development of the Model Equations 


In order to describe the non- conservative species transport of 
water borne constituents in Mobile Bay, knowledge of the current pat- 
tern and mixing characteristics must first be available. This informa- 
tion was developed in the study by Hill and April (12) titled " A Hydro- 
dynamic and Salinity Model for Mobile Bay” and is used to input veloci- 
ty and dispersion coefficient data for use by the NCSTM model. With 
this as background, the remaining portions of the chapter will be used 
to develop the NCSTM for Mobile Bay as applied to total coliform bacte- 
ria, BOD (biochemical oxygen demand^ and DO (dissolved oxygen). 
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The differential equations used in the Non- conservative 
Species Transport Model originate from the application of the law 
of conservation of mass over a differential element in space through 
which the liquid under consideration is flowing. Because of the 
shallow nature of Mobile Bay and relatively good mixing charac- 
teristics resulting from the interaction of fresh river water with 
seawater from the Gulf, the general equation can be modified to a 
two dimensional non-steady-state form. This equation can be adapted 
to describe the transport and fate of various non- conservative 
species by application of specific source and sink terms occurring 
at the boundaries of the system. In this study these models are 
referred to as the Total Coliform Bacteria, BOD, and DO models. 

3.2.1 Assumptions and Restrictions 

In order to derive an equation that will accurately predict 
bay system behavior while remaining solvable, a series of assump- 
tions and restrictions applicable to Mobile Bay are defined. These 
assumptions and restrictions are summarized in the following para- 
graphs . 

(a) Two Dimensional System 

As has been described in Section 3«1» the depth of Mobile 
Bay is very small (average 9*81 ft.) as compared to its length 
(approximately 31 miles) and width (ranging from 8 to 24 miles). 
Because of the effect of prevailing tidal action, the bay system 
as a whole can be considered vertically well mixed. Values of 
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the system variables at any point within the "bay can thus be 
considered a constant average value at any depth. The system 
can be reduced to a two dimensional one in which only changes in 
the longitudinal and latitudinal directions will be studied. 

(b) Tidal Cycle Average 

Data available for total coliform for verification of the 
NCSTM are collected on a spot sampling basis and do not represent 
within-tidal sampling. Because of this sampling method, verifi- 
cation of the model must conform to this pattern, i.e, a tidal- 
average basis. In all cases where current and dispersion coef- 
ficient are used by the NCSTM, tidal average values are computed. 
These values are subsequently used to calculate coliform distri- 
bution patterns representative of the data available for verifi- 
cation. Furthermore, these data are combined to form monthly 
average coliform concentrations to permit the analysis of the 
computed results. The NCSTM can be exercised on a within-tidal 
cycle basis provided that suitable data become available to 
permit calibration and verification on that basis. 

As a result of this method of solving the equation, the 
NCSTM becomes a quasi-steady state solution of the equation of 
change . 

(c) Constant Density and Viscosity 

Because of the interaction between seawater and fresh water 
in estuaries, density variations can exist. These density varia- 
tions are observed as salt wedges, bores and other phenomena 
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■which result in sharp discontinuities within the water masses. 
Mobile Bay is no exception to this rule; a salt wedge forms near 
the Main Pass and extends to various levels depending on the 
seasons and fresh water discharge rates. However, when consi- 
dering bulk fluid transport, density induced current and mass 
transport effects are normally negligible. Furthermore, when 
variations are averaged over a period exceeding the tidal cycle, 
they can likewise be neglected with little error introduced. 

Such is the case in this study. Monthly average mass transport 
distribution patterns are projected for total coliform, BOD and 
DO, In this model, density variations are considered negligible 
and are omitted from the model equations. Similarly, viscosity 
changes are also considered negligible, and the Newtonian law 
of fluid motion applies. 

(d) Binary Mixing and Variable Dispersion Coefficients 

The NCSIM considers species transport to be governed by 
Fick's Law. This is to say that the species in question forms 
one component while the rest of the water phase (including all 
other species) forms the second component of the system. There 
is no evidence indicating the effect that other water borne 
components have on mass transport of the components under study. 
In the absence of such information the assumption that the system 
behaves as a binary mixture will be adopted. 

The dispersion coefficients (E^Ey) in this study are 
affected by three elements, i.e. the turbulence of the water 
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column, the vertical mixing, and the tidal -cycle-averaging com- 
putation procedure. Because laboratory estimations of these 
coefficients differ greatly from field observations, the confir- 
mation of a set of dispersion coefficients that describe the 
mixing behavior of a system is difficult to obtain. In this 
study, the empirical equation developed by Holley, et al fo^) i s 
employed, which states that the dispersion coefficient is a func- 
tion of the bottom friction, the maximum current velocity over 
the tidal cycle, and the water depth. For Mobile Bay the bottom 
friction and depth are nearly constant; thus the change in current 
velocity outweighs the influence of the others, and becomes the 
controllig factor, The dispersion coefficient is therefore cal- 
culated using the maximum localized velocity over the tidal cycle. 
This correlation will be detailed later in this chapter. 

(e) Homogeneous Water Temperature 

In this study the water temperature of Mobile Bay is assumed 
to be constant at a unique temperature all over the bay. This 
assumption may introduce some inaccuracy for some locations 
within the bay where localized temperature gradients exist, e.g. 
the lower portion of the bay where seawater at a slightly diffe- 
rent temperature intrudes and causes moderate temperature inhomo- 
geneity. However, for a well-mixed tidal -smoothed model applied 
to Mobile Bay, the errors caused are negligible and the assump- 
tion of homogeneous temperature is reasonable. 


Based on the above restrictions and assumptions, the general 
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species continuity equation will be simplified to a form applicable 
for use in describing material transport in Mobile Bay. 

3*2.2 The General Species Continuity Equation 

Consider a differential element having length, width and height 
of x, y, z, respectively, fixed in space. Next consider the flow 
of a binary liquid into this volume containing species A with a con- 
centration of 



A differential 
element 


z 




The law of conservation of mass for this system, simply stated, is 


( rate of mass of A in ) - ( rate of mass of A out ) 

+ (rate of production of A by chemical or biological 
reaction or other sources other than by convective 
flow or diffusion ) 

= ( time rate of change of mass A in the element ). 

Therefore, the following quantities may be formulated* 

Input of A across face at x *. (NaxIx) A z 
Output of A across face at x-fAx s (Na x |x+ax) 

Input of A across face at y s (N^|y) axaz 
O utput of A across face at y+ Ay • (N^y jy+ Ay) Ax 
Input of A across face at z * (% K | z ) A x Ay 
Output of A across face at z+ az s (Na z | z + a2 ) ax A y 
where Na = mass flux. 

Rate of production of A by chemical reaction ( or any 
other generation and/or dissipation mechanism other than 
the advective flux term ) j r^ ax ajt az, 


Time rate of change of mass of A in volume element : 
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Substituting the above terms in the general mass balance equation, 
dividing by the differential volume ax Ay a z, and taking limits as 
ax, Ay, az approach zero, gives Eq. (3-l). 

dpA , ( dNAx , , 3^Az \ „ r (3-1 ) 

dt Vbx dy 32 / A 


The quantities N Ax , N A y, and N Az are the rectangular components of 

(?) 

the mass flux vector defined by Bird, Stewart and Lightfoot v 7 asi 

Na = <°A V A . (3-2) 

In vector notation, Eq. (3-1 ) becomes 

+(y’N A )= r A . (3-ia) 

Erom Pick’s first law. of binary diffusion^), 

% = w A ( N A t l B ) - ^ Dab VW A , (3-3) 

p 

where W A = ■■■■ * ~ mass fraction of A, 

r 

Dab “ mass diffusivity in the binary system. 

Equation (3-4) is obtained by substituting N A in Eq. (3-2) and 
transposing terms, 

— f- y*( W a ( ^ a + Nb )) ~ ^* ( ^ -Pa 6 ^ Wa) + ^ .... (3-4) 

where W A = (3-5) 

r 

% = Pa v A , (3-6) 

% = Pb (3-7), 

and V = mass average velocity =-~( f A V A + Vb)..(3“8) 
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Using Eq, (3-8) in expanded form, 

W A ( Ba+%) - ( ft 7 a + P B 57b ) (3-9) 

" (a ' (~p ( Pa Vf, + <°b Ve, )] 

= P A -V. (3-10) 

Eq. (3-4) can be rewritten as 

ff- + v ' (ftt 7 ) = v '(P^e viv A ) + r A (3-11 ) 

Expanding the divergence on the left hand side of Eq. (3-11 ) gives i 

+ft(v(7)4-(v-vf5) = v.(p D( , B viv 4 ) + r / , (3-12) 

If P , the overall density of the liquid system, is constant, then 
( v ■ v) " 0 and 

V .(f Dab v¥ A ) = V • (D AB v pw A ) (3-13) 

= v-(DabvPa), (3-14) 

and Eq. (3-11 ) "becomes 

~2t~ + (y- v ) = (y * d ar v ) + r A (3-15) 


This equation, expanded in rectangular coordinates, is 


(3-16) 



In this equation the instantaneous fluctuation in velocities 
and concentration with time (the turbulent phenomena) are not con- 
sidered yet. In estuarine flow where tidal action is a controlling 
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influence, turbulent effects are important. It is convenient to 


describe a turbulent variable by a time -smoothed 
tional deviation term, as illustrated for V 

term and a fluctua 
and Dae below j 

fA = ?A + (a* 


.. (3-17) 

V = V + V* 


■ • (3-18) 

- Bab + d ab ’ 


. (3-19) 


where barred variables are time smoothed parts and primed 
variables are fluctuational deviation parts. In Eq. (3-18), 
for example, 



= time smoothed V, ....... (3-20) 

where t 0 is a time interval which is large with respect 
to the time of turbulent fluctuation. 

Figure 3*2 shows this relation for the velocity Vj this figure can 
be equally applied to and d ab- If we take the time average of 
Eq. (3-16) by integrating each individual terms over the time interval 
t 0 and then dividing by t 0 , then all the fluctuational deviation terms 
Pa*» V and D^g' will vanish under integration. However, a quantity 
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such as V x * , the time-average of the product of two fluctuational 
terms, would not vanish under integration. In fact, it contributes 
appreciably to mass transfer, and can be represented with a diffusi- 
vity term which will be discussed later in this section. Attention 
will now be turned to the time -averaging of Eq. (3“l6). 

Assuming the term r^ in the right hand side of Eq. (3-16) can 
be represented by a first order reaction, we have 

+ -w tw + ( W + ( e A v 2 ) 

-hf A ( 3 - 21 ) 

where 

and overbars denoting vectors are dropped for the sake of simplicity# 
When V, are each replaced with Equations (3 -l6)f (3 _ 1?) an -d 
(3-18), we obtain after time-averaging 

-#-^^r(?/.V,) + -^r(P A V y )+^-(f A \4) 

- kif fl ( 3 - 22 ) 

Eq. (3-22) is the time-averaged species continuity equation in which 
V x * is the so called turbulent mass flux term. Over bars in 
Eq. (3-22) denote time -aver ages. 

By analogy with Fick’s law of diffusion 
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fA’ V - - D AB tx # 0-23) 

where is the "turbulent diffusivity". 

Substituting Eq. (3-23) in Eq. (3-22) and rearranging gives 

-$- + -k(W*) + ir&vy)+ 

= wi < Dab, + D, 6t , ) if] + £ [ ( B A6y + D«., ) -ff ] 


+ TF f < D m , + D*^) -ff-J - L ft (3-24) 

where an ^ ^ABtx ' be com " b l neci to gl ve a single term 

e x which is called the "eddy diffusivity”. 

<BlUS e x “ ®AB X + ^ABtx (3-25) 

e y “ -^ABy + ^AB-ty ...... (3-26) 

e z = %B Z + %B tz (3-2?) 


Combining Equations (3-24) to (3-2?) gives 


( f A v x ) - jy- c e A v, ) - jl- c v 2 > + a fc ) 


at 


+ + ir^ If) 


(3-28) 


Breaking the mass flux term, dropping the bars, rearranging, and 
loting that ( v»V ) = 0 for an incompressible fluid, Equation (3“28) 


jecomes 


at 


(Vx a x & + V; 


JL& + w -2 % 
ay v? az 


*—) + (& 


e x 


ax 


^ c 3 ft , 

3v + 


_afi_ \ _ 


Lr P. 


(3-29) 
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Eq. (3-29) is the three-dimensional general equation for the 
non -conservative species. However, in its present form, this equation 
is difficult to apply and requires numerical solutions which are 
lengthy and complex. Moreover, even if a solution to this equation 
could be obtained, the accuracy required for the initial conditions 
and boundary conditions would demand field measurements in excess of 
the capability of conventional field surveys. To circumvent this, the 
equation can be simplified to a two-dimensional form, then averaged 
over a tidal cycle to match existing field data for calibration and 
verification purposes. 

3.2.3 Simplification into a Two-dimensional Form 

According to the description cited in Section 3*2.1, the 
general equation for the NCSTM can be reduced to a two-dimensional 
form. This is done by neglecting the third-dimensional component of 
each variable and then vertically integrating Eq. (3-29) from the 
bottom to the surface of the bay water, then dividing the integral by 
the depth of the water column. For example, the depth-smoothed 
current velocity may be written as 

7x = ~tJ z b V x dz > (3-30) 

where - z at the bottom of the bay water, 
z s = z at the surface of the bay water, 

D = z g - = depth of the bay water, 

V x = x- component of the current velocity. 

Similar expressions can be written for Vy and other variables. 
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Eq. (3" 30) is exactly similar to Eq. (3-20 ) for the time- 
smoothed turbulent velocity, and a figure exactly similar to Fig. 3.2 
may be drawn for V x in Eq. (3"30) except that the horizontal axis must 
now be replaced with the depth of the water column. Another three 
equations similar to Equations (3-17) to (3-19) can be written for 
each of the variables V and e, and they can in turn be substituted 
in Eq. (3-29) » By doing so another set of non-vanishing Pa/V* terms 
will occur which, as in the time -smoothing operation, can be replaced 
by another analogy to Fick's law of diffusion. This gives rise to 
other diffusional mass flux terms which can again be combined with 
the eddy diffusivity terms to form a new set of diffusivity terms 
e x and e-y which are the so called dispersion coefficients. 

By assuming negligible variation in depth over the bay and by 
going through the averaging procedure similar to that from Eq. (3-22) 
to Eq. (3-29) » a two-dimensional species continuity equation is 
obtained: 

2 1 (V* 3* 3X ; 

+ ~ fr ( §y ~ , (3-31) 

where V x , and Vy are each vertically-averaged variables, and 
e x and §y now include the diffusional effect of vertically averaging 
the general species continuity equation. To suit the presently 
available data for the NCSTM, Eq. (3-31 ) must be further simplified 


into a tidal-smoothed form. 
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3.2.4 Tidal-smoothed Non- conservative Species Continuity Equation 


Ea. (3-31) is further simplified, "by averaging over the tidal 
cycle period T 0 , This is done "by a procedure similar to that in 
Section 3*2.2 for the time-smoothing of turbulent variables except 
that the time interval for integration now is the tidal cycle period 
T 0 which is much larger than t 0 . The variable V x , for instance, in 
smoothed form becomes 


Vx 


i_ ( t+T ° 

ToJt 


v x at 


(3-32) 


Similar expressions can be written for Vy and other variables. 


Here again it is noted that Eq. (3- 32) bears an exact 
resemblance to Eq. (3-20), and again a figure identical to Fig. 3*2 
can be drawn for V x , except that t is now replaced by T 0 (approximate- 
ly 25 hours). As was done in Section 3-2. 3» another set of expressions 
similar to Equations (3-17) to (3-19) can be written for the variables, 
and again a set of non-vanishing mass flux terms would occur,. They can 
similarly be replaced by an analogy to Fick's law of diffusivity. The 
newly created diffusivity terms can be combined with e x and ey in 
Eq. (3" 31) to form a new set of diffusional terms E x and Ey. Therfore 
by going through steps similar to those used in obtaining Equations 
(3“22) to (3-29 )» a tidal -smoothed two-dimensional non- conservative 
species continuity equation may be obtained! 




- ( V, + v y -If-) + £(E x -*&- ) 
+ -§y-(£y^)-kA , 


(3-33) 



35 


in which and Ey- are the ti&al-averaged dispersion coefficients. 
They contain the diffusional and dispersional mass transfer effects 
from time-smoothing, depth-smoothing, and tidal-smoothing the non- 
conservative species continuity equation. 

For the sake of rotational convenience, Eq. (3~ 33 ) can "be 
written into a more general forms 



-(u 


3X 


+ V 


3 ?a 
i>y 








+ (3-34) 


where U = net x- component current velocity over the tidal cycle, 

V - net y- component current velocity over the tidal cycle, 

E^r = x-component dispersion coefficient, 

Ey = y-component dispersion coefficient, 

2Sf = all the sources and sinks of the non- conservative 
species A, 

Equation (3 _ 3^) is the equation used in the NCSTM for the 
Mobile Bay system. This equation can assume different forms according 
to the difference in the term ZSj_, the mechanisms of generation and/or 
dissipation of the specific non -conservative species under study. 

3*3 Model Equations for Different Non- conservative Species 

Eq. (3-34) is applied to various non- conservative species, each 
having a distinctive mechanism of replenishment or consumption in the 
real system. This results in the total eoliform, BOD (biochemical 
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oxygen demand.) and DO (dissloved oxygen) models for Mobile Bay. 
3*3*1 Total Coliform Bacteria 


Total coliform bacteria group mean MPN has long been used by 
control agencies as an indication of the pathogenic bacteria content 
in waterways (-^6, 28) an ^ as a criterion for the certification of waters 
for the harvesting of shellfish, Ihe fecal coliform group, which is 
an indication of pathogenic bacteria derived from the excreta of human 
and other warm-blooded animals, has been recommended as a substitute 
standard for the certification of shellfish growing waters, and total 
coliform-fecal coliform relationship has been studied(^) , This 
relationship, usually expressed in the form of coliform-fecal coliform 
ratios, are subject to variations in the various bacteriological 
sources. Moreover, since all types of coliform organisms (fecal, 
non-fecal and intermediate) are found in feces, the absence of fecal 
coliform alone in waters designated for human use and contact is not 
a satisfactory criterion of acceptability. For the sake of safety, 
the standard test of the sewage pollution remains in terms of the 
total coliform group, although it has been argued that it is too 
stringent^) . In this study, the use of mathematical modeling as a 
predictive tool in determining the distribution of the total coliform 
group is studied for Mobile Bay. 

In studies of streams, the generation and dissipation terms 
for total coliform may contain: 

(l) Upstream runoff, 
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(2) Replenishment along the stream, 

(3) Reaction dissipation (die-off). 

In this study, the runoff term is expressed in terms of loadings 
(boundaries) at the mouths of rivers flowing into the bay, the repleni- 
shment term is neglected, and the die-off dissipation term becomes the 
main sink of total coliform bacteria. 

The coliform bacteria transported into waterways are investiga- 
ted and assumed to diminish by dying off at a rate proportional to the 
residual concentration, which is the same as a first order reaction for 
the stabilization of organic matter, radioactive decay, and many other 
natural phenomena . In equation form this is expressed ass 


-Jjf- = VB (3-33) 

The 2Sj_ term in Eq. (3-3*0 is thus 

r = - -|f- = - V B (3-36) 

where B - total coliform concentration in MPN/l00ml f 
t - time in days or seconds, 

Kp = dieoff rate constant in day“^ or sec~^. 

Substituting r in Eq. (3-3*0 gives 

. = — / u + w _?B . \ -f ( JL_ p j3B_ , ) 

v. u ax + v ay ) ' \ bx dx + zy By ' 

- (3-3?) 

Attention will now be turned to the correlations of dispersion coeffi- 
cients ^ and Ey and the dieoff rate constant K r . 
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Correlation for Dispersion Coefficient 

The dispersion coefficients E* and E y have "been studied by 

many workers, and several correlating equations have been derived 

through analytical treatment followed by experimental verification 
13) 

* Experimental results in natural streams, however, have not been 

within the expected range j deviations as large as several orders of 

magnitude have been found. In this study, the correlation by Holley, 
(13) 

et al . is adopted for the calculation of dispersion coefficients t 

S = 10 ° n s max R 5/6 . ....... (3-38) 

where E = dispersion coefficient, 

n = Manning's coefficient of bottom friction, 

S max « maximum absolute velocity over the tidal cycle, 

R = hydraulic radius 

_ cross sectional area of f low 
wetted perimeter 

When R is in ft. and Sj^ is in ft. /sec., E is in ft 2 /sec. 

In the case of Mobile Bay^ 12 ^, 

n - O.Olf to 0.018, 

R * 0.5 D, 

where D = average depth of the bay. 

From Eq. ( 3 - 38 ), 

E# 4.024 to 4.080 times S Ina:>c . ....... (3—39 ) 

The dispersion coefficient is therefore a linear function of the 
amplitude of the tidal velocity. The E value best suitable for a 
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certain species is then obtained by multiplying maximum tidal veloci- 
ty calculated by the Hydrodynamic Model by a correction factor const- 
ant (C.F,), that is 

E = (C.F.) • (3-40) 

In the calibration period, this factor is found by interfacing the 
dispersion coefficient with other model parameters; various values are 
used until the calculated results match the actual data. Different 
species behave differently' in the physical system, therefore it is 
expected that different values may be required for the calculation 
of other substances. 

Correlation for the Dieoff Rate Constant 

The dieoff rate constant, K r , as in any first order chemical 
reaction, is a function of temperature. Surveys performed by many 
workers on a number of rivers have given a range of K r values from 
0.26 to 0.46 in cool weather, and a range from 0.46 to O .96 in warm 
weather^). These ranges are adopted in this study, since field 
data to establish the reaction rate coefficients for Mobile Bay are 
non-existent. 

The temperature dependence of K r is expressed in the form 
k t = k 2O -0 t - 20 ... (3-4i) 

where Kq» = K at any temperature T°C in day“^ or sec"-*-, 

K 20 “ K at 20°C in day""-*- or sec"-!', 
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0 = a constant characteristic of the reaction, 
dimensionless, 

T = temperature °G. 

As an approximation, O .96 day - " 1- is designated as corresponding 
to 86 °F or approximately 30°C, and 0.26 day - '*' is designated as corres- 
ponding to 50°F, which is 10°C» These two K's give the value of 0 as 
1.06? and the value of Kgo as 0,50 day -1 . Thus 

K t = 0.50 x (l.06?) T-2 ° (3-42) 

is used for calculating th'e death rate of col if or m in this study. 

Assumptions and restrictions specific to this section are 
summarized below: 

(1) A first-order dieoff mechanism is assumed, for the coliform 
bacteria. 

( 2 ) The reference temperature for the dieoff rate constant obtained 
from literature is arbitrarily assumed to apply (due to lack of 
information) . 

(3) Replenishments along the perimeter of the Bay are ignored. 

(4) Upstream runoff is considered in terms of river discharges 
associated with a certain pollutant concentration, and in terms 
of boundary conditions at any apparent loadings. 

3*3*2 Biochemical Oxygen Demand 

BOD (biochemical oxygen demand) is the amount of oxygen 
required by bacteria while stablizing decomposable organic matter 
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under aerobic conditions. The decomposable organic matter can serve 
as food for the bacteria, and energy is derived from its oxidation. 
The BOD test is widely used to determine the pollutional strength 
of sewages and industrial wastes in terms of the oxygen that they 
will require if discharged into natural waterways in which aerobic 
conditions exist. 

Studies of the kinetics of BOD reactions have established 
that, like the dieoff of coliform, the reactions are first order in 
character (30), i. e . 

the rate of the reaction is proportional to the 
amount of oxidizable organic matter remaining at any time. A second 
order reaction mechanism has been under study for systems having 
critical oxygen deficiency. Positive results have been reported 
in the literature ( 32 ), 

Nevertheless, under the present bay condi- 
tions, the assumption of a first order reaction mechanism has been 
confirmed by most of the studies in other similar systems and is 
regarded as the standard practice. 

Carbonaceous BOD and Nitrogeneous BOD 

Figure 3-3 shows a typical BOD, or oxygen use curve, which is 
typical of laboratory BOD tests. This figure shows that there are 
two stages of BOD reaction, i.e. the Carbonaceous BOD and the Nitro- 
geneous BOD. Extensive studies have shown that the bacteria derived 
from soil or domestic sewage are actually a mixed culture of orga- 
nisms corresponding to large numbers of saprophytic bacteria (and 
other organisms that utilize the carbonaceous organic matter) with 





a certain amount of autotrophic bacteria, particularly nitrifying 
bacteria, which are capable of oxidizing noncarbonaceous organic 
matter (^3), phe nitrifying bacteria are found to be usually present 
in relatively small amounts in untreated domestic sewage. However, 
they are capable of reproduction*, their reproductive rate is small 
enough so that their population does not become sufficiently large 
to exert an appreciable demand for oxygen until 8 to 10 days have 
elapsed in regular BOD tests. In stream and estuary systems, their 
presence is affected by the nature of the waste material j field 
surveys are required to find the amount and rate of reaction for 
nitrogen eous BOD. Total omission of it as an important input can 

only be justified if a time of passage of less than 8 to 10 days 
at 20°G (or the equivalent time period at other temperature) exists. 
Therefore it is always safer to assume the coexistence of NBOD 
(nitrogeneous BOD) and CBOD (carbonaceous BOD) in'any system with 
complex flushing characteristics. 

The generation and dissipation terms for BOD, that is the 
” 2Si M term in Eq, (3”3^)» should contain 

(a) Replenishment along the watercourse (source), 

(b) Input from upstream runoff (source), 

(c) Resuspension from the benthic layer (source), 

(d) Deposition or sedimentation into the benthic layer (sink) , 

(e) Oxidation reaction use (sink). 

Eq. (3-34) may now be written for CBOD as 
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and. for NBOD as 
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( Kn + K 3n ) N + N Rl 



dN v 

^y ' 


(3-44) 


where L, N = concentration of CBOD and NBOD respectively in 
mg/liter, 

K^, = reaction rate constant of CBOD and NBOD 

respectively in day”* 1 ' or sec"\ 

K^» K-^ n = rate constant of resuspension and sedimen- 

-1 -7 

tation in day or sec x , 

Lp, Nr = replenishment along the watercourse in 
mg/liter-day. 


and K^ n are sometimes grouped into K m ; and are sometimes 
grouped into K r ^. Values of K-^ = 0.34 day *** and K-^ — 0.14 day 
Both at 20°C, are used by the Galveston Bay study (5> 6). Again, as 
in the case of total coliform bacteria, the replenishment terms Lr 
and Nr are assumed negligible for both CBOD and NBOD. Upstream 
runoff is also expressed in terms of loadings at the mouths of rivers 
draining into the bay and appears in the boundary conditions. 


It would be ideal if there are enough informations on all 
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the entities involved in the above equations for Mobile Bay or for 
other similar systems. When niether is available, specification is 
done on the basis of similar behavior in waterways and streams. 

In some stream studies^), values of ranging from 0.49 to 
3.5 day are used, and values of ranging from 0.1 to 2.5 day 
are used. At 28° C, /k^ ratios of 2.375 and 2.362 have been used 

for streams flowing at 0.922 and 0.510 ft/sec respectively; values 

—1 “1 

of O .76 to 0.95 day” x and K n values of 1.9 to 2.5 day” were used at 

28° G. These values are extrapolated to 20° C with the expression 

k t = k 20 e T " 20 or k 2 o = k t e 20 " T (3-^5) 

for the BOD reaction, 0 =1.03 is suggested (6), 

Values of and 

K n , both at 20°C, are correlated as 


Kd,20 88 0*68 ± 0.08 day" 1 (3-46) 

and K n> 20 = 1*74 ± 0.24 day" 1 (3-47) 


These values can be adopted in place of K r< i and respectively, as 

a first approximation. Eq. (3-45) can then be used to extrapolate 
to temperatures other than 20°G, In adoption of information from 
other systems, it is assumed that the aquatic ecosystems from which 
information is solicited behave similarly to those in the Mobile Bay 
system. This is an approximation, and a trend analysis of the model 
results can be made. 

3*3*3 Dissolved Oxygen 


DO (dissolved oxygen) in waterways is important to aerobic 
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aquatic lives as atmospheric oxygen is important to men. Severe 
deficiencies of DO in water often result in fish kills. Therefore 
it is required that DO levels be maintained to support aquatic lives 
in a healthy condition at all times. Most of the critical conditions 
related to DO occur during the summer months when temperatures are 
high, rates of biological oxidation increase, and DO contents 
decrease to minima. Fig. 3*^ shows a solubility curve for dissolved 
oxygen in water saturated with air at 1 atm. The saturation solu- 
bility of oxygen is usually used in estuarine DO studies. 


Sources and sinks of DO arei 

(1) Surface reaeration (source), 

(2) Photosynthesis generation (source), 

( 3 ) Upstream runoff (source), 

(4) Biochemical oxidation demand (sink), 

(5) Benthic layer uptake (sink), 

(6) Respiration use by all aquatic lives (sink). 
Eq« (3-34) thus becomes, for DO, 


6B 

St 


( u M + v S 
' Sx By > B x 


a <^Hr> + 


ay 




Ki L + K 2 (D s - D) + P - R - S b 


(3-48) 


where D - DO concentration in mg/liter, 

D s = Saturation solubility of oxygen in water in 
mg/liter, 

Kq = rate, constant of biochemical oxidation demand 
in day“^, 
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Kg = rate constant of surface reaeration in day - ^-, 

L = BOD concentration in mg/liter , 

P - Photosynthesis generation in mg/liter-day, 

R = Respiration use in mg/liter-day, 

= uptake by benthic organisms. 

The following paragraphs will discuss various terms in Eq, (3~^8). 
Biochemical Oxidation Sink K-] L 

Ibis sink is the consumption of DO by BOD. Strictly it should 
be written as K^L + K^N instead of KqL to account for CBQD and NBOD 
separately. 

Surface Reaeration Source 

This source and the photosynthesis generation are considered 
the primary sources for DO. The reaeration is regarded as first 
order as it is in many similar studies. At atmospheric pressure the 
rate constant Kg is considered as a function of temperature and other 
physical effects. In laboratory studies for obtaining Kg, the 
temperature effect is first fixed by finding Kg at 20°G. Values of 
Kg at other temperatures are then extrapolated using an equation of 
the form similar to Equations (3-4l) and (3-45). 

Studies by many workers on various streams and estuaries have 
resulted In the following empirical equations for Kg^^). 
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K 2 = 


5.026 v 0,969 
H l76?3 


(3-49) 


k 2 = 3.739 -~^rry 


K 2 - 0.00125 ( 1 + N f °*3 )JI 


gs_ 

H 


K. 




H' 


T^y 


k 2 - 


( % v ) 0,5 


2.303 H 


T3 


(3-50) 

(3-51) 


(3-52) 

(3-53) 


Of these, Eq. (3—53) by O'Connor and Dobbins is by far the one most 
often used. It is recommended for use in this study "because of its 
consistency in dimensions and covenience in use, In Eq. (3“53)i 

V = stream velocity, 

H = depth, 

Djvj = molecular diffusivity of oxygen 
= 0.81 x 10"** ft 2 /hr at 20°C. 

Kg at any other temperature T is calculated by 

K 2,T°G = k 2,20°C © T ~ 20 » (3-54) 

where a suggested value of 0 for the DO reaction is 1.02^\ and 
T is in °C. 

Photosynthesis Generation P and Respiration Sink R 


Recent studies (20 ) 


show that oxygen contribution by photosyn- 
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the tic activity P is a primary source of DOj its value can predomi- 
nate , or be equal to the respiration sink R, or be smaller than R 
under different circumstances, and should be experimentally evaluated 
instead of being stochastically neglected by assuming an gross equi- 
valence with R. 

The photosynthetic rate P is a function of radiation intensity 
and the phytoplankton population, which can in turn be functions of 
time, temperature, depth, and position, 

A time-varying P of the form of a half -cycle sine wave is 
suggested by 0* Connor and Di Toro^^s 

P(t) ~ P m sin [~_(t-t s ) ] when t s ^t^t s + P .... (3-55) 
= 0 when t s +P<t^t s + 1 .. (3-56) 

where P(t) = time varying rate of photo synthetic oxygen 
production in mg/liter-day, 

P m - maximum value of P(t), 

t g ~ the time at which generation begins in days, 

P = the fraction of the day over which photosynthesis 
exists . 

The periodic expression of Eq. (3-55) can "be expressed as a 
Fourier series and used for the long-term effect of photo synthetic 
oxygen generation, or can be used in the time-varying DO model to 
calculate DO at different times within a tidal cycle* It can also 
be integrated over a tidal cycle to provide an average term P av , as 
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is of interest in this study due to the lack of time- varying field 
data. The respiration term, unlike the photosynthesis term, is 
assumed to be constant over a certain period. The temperature 
effect on respiration can be expressed in the equation 

R T = V e rT (3-57) 

where Rp = respiratory rate at some temperature T, 

Rq = respiratory rate at 0 °C, 
r = constant to be determined by experiment, 
e = base of the natural logarithm. 

Studies made by Riley (9) on Long Island Sound found R 0 for 
winter and for summer to be 0.020 and 0.015 mg. of carbon consumed/day/ 
mg. of phytoplankton carbon respectively. Conversion of units is 
required in adopting these values in the DO model. 

Wright tabulated monthly averages of P and R for various 
streams during different months (April to October, 1957 and April to 
October, 1958) as a function of phytoplankton densities. These values 
can be adopted before more suitable data become available for Mobile 
Bay. 


Benthic Uptake 

Although in streams this term is often neglected by assuming 
bottom scour due to high speed of flow, this term deserves more 
consideration in an estuary like Mobile Bay. However, data on this 
sink are not available. In the Galveston Bay Study the equation 



( 3 - 58 ) 
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S b = 2.0 e°-07(T-20) 


day 


was adopted due to lack of data for Galveston Bay as of the time of 
report issuance. In this study, the "benthic uptake can "be (a) negle- 
cted by assuming that is is mainly due to benthic bacterial respirati- 
on (21) and Is included in the respiration, term R (the nature of R data 
adopted have to justify this), or (b) calculated using Eq. (3 "58 ) with 
conversion of gm/m^.day into mg/liter -day by incorporating the local 
depth of bay cells. 

With the development, simplification and adaptation of the 
general non- conservative species continuity equation to total coliform, 
BOD, and DO completed, attention is now turned to the numerical method 
used to effect solutions. This will be followed in the next chapter 
by a discussion of the calibration and verification methods used to 
test the coliform model. Results of the coliform model including 
parametric studies involving varying river discharge rates, wind 
conditions and temperatures are presented and discussed in Chapter V, 

3.4 Numerical Solution of the Non-conservative Species Equation 

Finite difference equations can be written for the partial- 
differential equations developed for the various models in the 
preceding sections. Finite space increments, ax and Ay, and a 
finite time increment At, are selected based on the stability criteria 
insuring a correct solution. As shown in Figure 3 *5* a grid system 
consisting of 38 rows and 21 columns formulating 798 square grid 
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Fig. 3.5 The grid system (21 x 38) superimposed on the Motile 
Bay system for the finite differencing technique. 
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cells of 2 km. by 2 km., is superimposed on the Mobile system. The 
computer program for the Salinity Model developed by Hill and April 0-2) 
is adopted and modified to implement the finite difference equations 
for the NCSTM for Mobile Bay. Initial and boundary conditions are 
supplied and grid to grid computations are effected until relaxation 
occur, that is, until results computed for two consecutive sweeps of 
the grid system are within tolerable differences. 

3,4.1 Finite Difference Techniques 


Forward, backward, and central differences have been used in 
finite differencing methods for solving ordinary and partial differen- 
tial equations. Their basic forms are summarized below. 

Upon proper selection of ax and Ay W, 


forward difference 


M. „ x {l +i) - x(i ) (3 _ 59 ) 

dy Ay 


backward difference 


dx x(l) - x(l-l) 

r — — - 

dy Ay 


(3-60) 


central difference 


(3-61) 

dy 2 Ay 


By properly subscripting the dependent variables, the finite diffe- 
rence equations can be written to include values from both neigh- 
boring columns and neighboring rows, that is, the two dimensions 
x and y. In this study, the forward difference is used for first 
order derivatives, and the combination of forward and central diffe- 
rence formulas is used for second order derivatives. 
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3.4.2 Finite Difference Equations for the NCSTM 


Application of the finite difference method is made here to 
the computation of total coliform concentration in Mobile Bay, 
beginning with Eq_. (3-37) rewritten here as Eq. ( 3 - 62 ). 


3B 

at 




SB 


dy 


- K r B 


(3-62) 


Upon finite differencing, 


jJB_ _ B 1 (i, J) - B(l,J) 

at ~ At 


3B_ „ B(I+1,J) - B(I,J) 
8x ~ AX 


a® ~ b(i, j+i) - b(i,j) 

dy ~ Ay 


U(l+l,j) - U(I,J) 


V(l,J+l) - V(I,J) 
2 


(3-63) 


(3-64) 


(3-65) 

(3-66) 

(3-67) 


where B(l,j) = B at time t, B’(l,j) = b at time t + At, 
in cell (l, J) . 

Furthermore, by the product rule, 


_d 

ax 

_a 

dy 
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dx ^ Bx ^ 
“ dy ^ By * 
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3 2 B ■> 
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3 2 B >. 
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( 3 - 68 ) 


(3-69) 
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For the x-component, 

SE „ E(I+1,J) - E(I,J) 


Sx 


AX 


(3-70) 


3B „ B(I+1,J) - B(I,J) B(I,J) - B(I-1,J) 

ax ~ AX ~ 2Tx *•' 

B(l-H,j) - B(l, j) B(I,J) - B(I-1,J) 


a 2 B 

Sx2 


ax 


Ax 


AX 


= B(I+I,J) - 2B(I,J) + B(I-1,J) _ ... (3 _ 72) 

Combining Eq. ( 3 - 68 ) with Equations (3-?0), (3-71) and (3-72) gives 


dE 

dx 


5B 2r f E(l+1, J) - E(l t j) 1 . r B(l+l f J) - B(I,J) 
C>X ~ l AX J ^ AX 


+ E(I,J) 


B(I+1,J) - 2B(l,j) +B(I-1,J) 
( ax) 2 


E(I+1,J) - S(I-1,J ) ] f B(l+ l, J) - B(I,J) 


2 AX 


A X 


+ E(I,J) 


B(I+1,J) - 2B(I,J) + B(l-1, j) 


(ax) 2 

2( 5c) 2 E(l+1 , J) [B(I+1,J) - B(I,J)] 

E(I-1,J) fBd+l.j) - B(I,J)] 

2( ax) 2 L J 


+ e(i.j)- 77^)2 - zbCaJ) +B(I-1,J)J 

— B— r E(I+1,J) [b(I+1,J) - 

2(ax) 2 
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- - ■ - . E(I-1,J) [b(I,J) - B(I-I.J)] 

2 ( 

+ — — E(I,J) [b(I+1,J) - 2B(I,J) + B(I-1,J)] 

(iXf 

- , 1 ,, E(I+1,J){[b(I+1,J) - B(I,J)] 

2( Ax)^ 

- E(I-I.J) [b(I,J) - B(I-1,J)} 

+ 2E(I,J) [B(I+1,J) - 2B(l, j) +B(I-1,J)]} ... ( 3 - 73 ) 
Similarly for the y-component 

17 ( Ey “If-) * ^7^2{ E(I ’ J+1) ( B(I ' J+1 ^ ‘ 

- E(I,J-1) [b(I,J) - B(l,J-l)j 
+ 2E(I,J) [b(I,J+ 1) - 2B(l, j) + B(l,J-l)jj 


( 3 - 7*0 
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Finally all these finite difference formulas are substituted into 
Equation ( 3 - 62 ) to give 


At 


--ft 


U (I + >,j) - If (1,1) 

z 


BCl-H, J) - 6 Cl 
AX 




+ yiLjjtji - 7(^J) j j" g(X,3 + i)- 8(^,3) jj 

+ { eo-xi.j) [b(i+i,j;- bcx,?;J 

- id-1,7) - 8(x-/,J)J 

+ 2 £(i t j) [sa-t i,j; - 2 B(i, i) + 6 (i-i, j;J| 

^ A { (l, J+/^ — & ( 2 , 

- E(i,mj - b<xj-oJ 

+ 2E(l,?)[e (1,7+1) -2B(I y 7) +8(2, 

~ krB(1 J 3). ( 3 - 75 ) 
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Rearranging and solving for B*(l,J) results in the desired equation 
to be applied to each grid cell in the bay, 

+ Twy. { E(I*UJ [B(I+i, 3 ) - B(x, 3 jJ 

- £(x-i, j) • - Ba-1,3;} + 2 e(i,j;[b(i+i,j; 

- 20d, j) + b 

+ E ( - B(x, J)) 

+ 2 £ (1,3.) • [b - Z 8(1,7) + 

- (At) KrRU.J) 


This finite difference equation is used to implement computer solution 
of the species continuity equation for total coliform in Mobile Bay. 
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3.4.3 Application to Mobile Bay 

Details of the development of the computer program have been 
cited in the work of Hill and April They not only include 
specific derivation of the species continuity equation for conservative 
species, hut the development of equations needed to specify the current 
distribution in the bay which is a critical input to this study. The 
specific aspects from that study which apply to the NCSTM are summariz- 
ed below. 

Finite Increments 


There are specific limitations on the sizes of the finite 
increments for the finite difference solution to be stable or to 
converge. For the species continuity equation^-^) these limitations 
are 


A X 




At 


< 


< 


<: 


^max 

2E^ 

Xmax 

( as ) 2 

2(Ex + %) 


(3-77) 

(3-78) 

(3-79) 


where as = ax = Ay in this study. 

The spatial increments in the x and y directions (ax and Ay) were 
chosen to be 2 km ( 6561.68 ft) each (see Fig. 3*5 on p. 33). A 
time increment At of 240 seconds was chosen to insure stability. 
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Boundary Conditions and Initial Conditions 

Boundary conditions (concentrations for certain border line 
cells) and initial conditions (concentrations for all cells at time 
t = 0) must be specified in order to solve the partial differential 
equation describing a system using a finite difference technique. 

For Mobile Bay all the land cell concentrations are set equal to zero 
and the partial derivatives at water-land connecting cells are set 
equal to zero. The concentrations of cells on the Gulf front of the 
grid system are set equal to zero in the coliform model. This is 
reasonable because the coliform bacteria levels become negligibly 
low in Gulf water. 

The Mobile River and the Tensaw River are each simulated with 
an idealized channel ten grids long flowing from the north to the 
south draining into the bay. Upstream runoffs from the rivers are 
expressed as boundary conditions at the mouths of the rivers where 
they flow into the bay. Values of these boundary conditions are 
adopted from data collected at water quality stations corresponding 
to or located near the boundary cells under consideration. Grid 
cells in the bay near suspected outfall of waste as reflected in 
field data are also assigned as boundary cells. Alabama Port (see 
Fig. 4.l) is an example of this behavior as reflected by the high 
coliform levels in waters adjacent to it. Therefore a boundary cell 
is assigned to the bay water near Alabama Port where the concentra- 
tion is fixed in performing the computation. The value of the 
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concentration is taken from a station corresponding to that cell. 

Initial conditions are set equal to zero for all the cells 
for the first computation when prior knowledge of the system beha- 
vior cannot be estimated. The first computed results are then stored 
into a data file. Provisions are made for subsequent computations 
to utilize the previous result as initial conditions. This serves 
to conserve some computing time which is important for this kind of 
calculation. 

Sources of Data 

In order to exercise the NCSTM, specific information from 
various sources must be supplied. These data and their sources are 
summarised below. Their formats are listed later in the appendices. 

(1) Monthly average river flow rates for Mobile River, the main 
source of pollution of Mobile Bay, are provided by the Alabama 
State Geological Survey at Tuscaloosa, Alabama, Data for the 
period January -August 1962 are selected for use in the coliform 
study to match the available coliform data. River flow rates 
for other periods starting from August 1928 are also available. 
The Geological Survey also provides daily river discharge rates 
for the Alabama River and Tombigbee River which can also be 
valuable to the short-term water quality study of Mobile Bay. 

(See Appendix B2,p. 181 ) 

(2) Wind data are obtained from the climatological data collected 

by the Weather Bureau of the Department of Commerce. Prevailing 
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wind speed and direction for each month are used for the study 
of total coliform. (See Appendix B3, p. 184) 

(3) Atmospheric temperature data over the Mobile Bay area are collec- 
ted by the U. S. Weather Bureau. Bimonthly average bay water 
temperature profiles of Mobile Bay have been compiled and pre- 
sented by Bault^). The latter forms the basis for the deter- 
mination of water temperatures to be used in the study of total 
coliform. (See Appendix B3, p. 184) 

(4) Total coliform data are provided by the Alabama State Department 
of Health for the period January -August 1962. Numbers of data 
points for each station range from two to five per month. They 
are averaged on a monthly basis to be utilized in the study on 
total coliform distribution in Mobile Bay. (See Appendix Bl, 


p. 169) 

Data specific to the hydrography of Mobile Bay are necessary 
for the NCSTMj they are adopted directly from the work by Hill and 
April Formats of input and output data of the NCSTM, computer 

program listings, and descriptions of the model variables are summa- 
rized in Appendices A1 to A4» 


General Computation Procedures 

The above data are used in the NCSTM to obtain the total 
coliform profiles for Mobile Bay in the subsequent chapters. Other 
water quality species can be investigated provided that pertinent 
field data are available. The general procedures of computation of 
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the NCSTM are summarized, "below: 

Step 1. River discharge rates, wind speed and direction, and other 
hydrographic data of Mobile Bay are fed to the Hydrodynamic 
Model of Mobile Bay * The Hydrodynamic Model calculates 
the net current velocities over one tidal cycle and the 
maximum velocities for each of the water cells of the bay. 
These data are then written into data files and stored in 
the memory of the computer. 

Step 2. Temperature, boundary conditions of total coliform concen- 
tration, and pertinent hydrographic data of Mobile Bay are 
fed to the NCSTM. Starting with zero initial concentrations, 
the NCSTM reads in the data file created by the Hydrodynamic 
Model and computes the total coliform concentrations for 
each cell of the bay. 

Step 3. Each pair of consecutive computations are compared until 
the computed results converge, that is, when the results 
yielded by two consecutive "sweeps” over the grid system 
are within acceptable deviation. In the computations per- 
formed in the following chapters, calculated results usually 
converge to within + 1 %. 

The final results are then compared to the field data for the 
purposes of calibration and verification. The following chapter deals 
with the verification of the NCSTM for total coliform bacteria dis- 


tributions in Mobile Bay. 



CHAPTER IV 


CALIBRATION AND VERIFICATION 

Total coliform group concentration data for various locations 
in Mobile Bay were collected by the Alabama State Department of 
Health for the period from January 1962 to August 1962. Figure 4.1 
shows the locations of the coliform sampling stations in Mobile Bay 
(i4). 

Figure 4.2 shows the grid locations with corresponding station 
numbers at which coliform concentration data are available. It is 
these data that are used for the purpose of verification of the 
Non -conservative Species Transport Model for Mobile Bay. These 
coliform group concentrations are obtained by analysis as described 
in the outline entitled "The Significance of EC- Positive Organisms 
in Gulf Shellfish Growing Waters" (see Appendix Bl, p. 170). 

The model is verified on a monthly basis, i.e. monthly ave- 
rage conditions are used, and the model results are tabulated and 
compared to the monthly average values of actual data. The 70% 
confidence ranges of the actual data are also tabulated to indicate 
the range in the monthly field data averages. The criterion for 
model verification is based on how well model-predicted results 
fall within the field data range at the several locations within 
the bay for any given monthly period. 

It will be shown in the following sections that the model 
predicts resonable results as compared to the measured data. 
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However, it should be indicated that collected field data used in 
the verification were not specifically obtained for mathematical 
modeling purposes, and thus represent selections based on availa- 
bility. A more detailed verification program, including synoptic 
data collection specifically in support of mathematical modeling 
efforts, would be required to confidently use model predicted re- 
sults for trend analyses on less than monthly frequencies. 

Additionally, there are no sources of data of any magnitude 
to suitably verify the biological oxygen demand (BOD) and dissolved 
oxygen (DO) models. As a result, the verification of the NCSTM is 
based solely on total coliform group concentration data. Extrapo- 
lation of conditions, using the experience gained from the total 
coliform model and available literature surveys concerning those 
concepts and laws governing the calculation of BOD and DO concen- 
trations in estuarine waters, can be used for preliminary trend 
studies. However, verification studies of the BOD and DO models 
will have to be made, including design and implementation of sui- 
table data collection programs in support of mathematical modeling, 
before the BOD and DO model results can be used in a truly predic- 
tive capacity. 

A.l Calibration and Verification Procedures 
Interaction with the Hydrodynamic Model 


Because of the dependence of the species continuity equation 
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on the hydrodynamic model of Mobile Bay for current distributions 
and dispersion coefficients, the first step in the verification 
procedure involves specification of data necessary for the proper 
description of the hydrodynamic behavior of the Bay. This includes 
the calculation of monthly average river flow rates, wind conditions 
and tidal conditions for the period for which total coliform group 
concentration data are available. 

Monthly average river flow rate data are available from the 
Alabama Geological Survey at Tuscaloosa, Alabama for the Mobile 
River near Mount Vernon, Alabama. These are reasonably split into 
two parts for the river discharges of the Mobile River and the Tensaw 
River which empty into Mobile Bay in the north. River discharge 
rates for Dog River are varied between 500 cfs and 5000 cfs, depen- 
ding on the month in which the model is to he exercised. Values 
used for verification studies during the period January to August 
1962 are shown in Table 4,1. 

Wind conditions, including speed and direction, are calculated 
as statistical averages for each monthly period during 1962 . These 
data are obtained from climatological data provided by the U. S. 
Weather Bureau (^7), Wind speed in knots and wind direction in de- 
grees from the x-axis are listed in Table 4.2 for the period January 
to August 1962 . 

The tidal cycle conditions are described by equations deve- 
loped by Hill and April for each of two locations where the bay 
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Table 4.1 River Discharge Rates for the Period 
January to August, 1962 in cfs 

Month Mobile River Dog River Tensaw River 

January 130,000 5,000 73,800 

February 100,000 3,000 51,700 

March 90,000 2,500 53,100 

April 130,000 4,000 56,900 

May 20,000 2,000 18,500 

June 15,400 1,500 10,000 

. . . ., — — — — — . - - — ... - .... 

July 10,000 1,000 9,200 


August 


8,000 


5 00 


4,500 
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Table 4.2 Temperatures, Dieoff Rate Constants, and Wind 

Conditions for the Period January to August, 1962 




Temperature 

Dieoff Rate 

Wind Conditions 

Month 

Constant K r 

Speed 

Direction 


°F 

day 1 

knots 

from 

0 deg. 

January 

49.5 

0.26 

12.3 

N 

90.0 

February 

53-2 

0.29 

12.0 

s 

270.0 

March 

61.3 

0.39 

12.6 

N 

90.0 

April 

67.9 

0.50 

10.7 

SSE 

292.5 

May 

78.1 

0.?2 

7.9 

SW 

225.5 

June 

81.4 

0.81 

5.7 

ME 

45.0 

July 

83.7 

0.88 

5.9 

SW 

225.5 

August 

84.2 

0.90 

5-2 

ENE 

22.5 
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interfaces with Gulf of Mexico waters, These equations describe 
the tidal level at the Cedar Point and Dauphin Island-Gulf boundaries 
and are represented as: 

HDI = 1,090 + 1.295 * cos (0.004188 • t + 0.0567114) 

HOP - 1.089 +1.177 * cos (0.004188 • t + 0.0032453) 

Unless there is evidence of conditions altering tidal behavior in 
the Gulf (i.e. storms, diurnal periods, etc), it is assumed that 
normal tidal conditions prevail over the monthly cycle. 

Using the above data as input to the hydrodynamic model, the 
corresponding output, including tidal cycle average velocities and 
dispersion coefficients for points within the bay, provides a des- 
cription of the period for which total coliform group data are avai- 
lable. 

Mon- conservative Species Model for Coliform 

Specification of inputs for the NCSTM for total coliform 
includes two types of information which are classified as follows : 

(1) cell data which includes velocity distributions and dispersion 
coefficients for each grid of the model as calculated in the 
Hydrodynamic Model; temperature data used in the calculation 
of the total coliform dieoff rate constants K r ; and 

(2) boundary and initial conditions of total coliform concentration 
data dealing with specific Inputs at the spatial and temporal 
limits of the model. These inputs are discussed below. 
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Cell Data and Conditions 

Net velocities over one tidal cycle for the grid cells are 
used as the x-component (u) and y- component (v) velocities in the 
Coliform Model. Maximum velocities over the tidal cycle for the 
grid cells are used to calculate dispersion coefficients (Ej, and 
Ey) according to Holley’s correlation equation (3-^0). The dis- 
persion coefficients thus calculated are used, after being modified 
by a correction factor suitably defined by monthly average field 
data during model calibration. These modified dispersion coeffi- 
cients are selected to provide the best description of the macrosco- 
pic mixing characteristics for the given species and conditions 
that exist within the bay. 

Additionally, the total coliform dieoff rate constant K r used 
in the model is calculated as a function of monthly average water 
temperature of the bay according to equation (3-^2). These tem- 
peratures are estimated from the bimonthly average water temperatures 
of Mobile Bay compiled by Bault^ 1 ). It is recognized that water 
temperatures are not uniform in the bay. The degree of mixing that 
occurs between sea water and river water within the bay will affect 
the temperature distribution. In this study temperatures are con- 
sidered homogeneous throughout the bay. Temperatures can be adjusted 
linearly between the values corresponding to Gulf of Mexico water 
temperature and river water temperature to approximate real system 
behavior. In this study where monthly average values are investi- 
gated the sea water intrusion effect can be neglected. This point 
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will be explained later in the discussion section (Section 4,3) of 
this chapter. 

Initial and Boundary Conditions 

Initial and boundary conditions for total coliform group 
concentrations are specified as described in Section 3-4.3 on page 6l. 

Cell Loadings and Dilution Factors 

Total coliform group concentration data for points recognized 
as having severe pollutant input into the bay are used as loading 
concentrations at each relevant grid cell as shown in Fig. 4.2. 

They are held constant throughout each computation. These points 
include (l) the mouth of the Mobile River leading into the bay, 

(2) the mouth of the Tensaw River leading into the bay, ( 3 ) the 
mouth of the Dog River leading into the bay, (4) the water adjacent 
to the Alabama Port, ( 5 ) Cedar Point, and (6) the mouth of the 
Bon Secour River leading into the bay. Loading at the Mobile River 
has been found to be the main source of pollution of Mobile Bay do). 
Values of the cell loading total coliform concentrations are shown 
in Table 4.3- 

Loading concentrations for the mouth of the Mobile River, 
when directly taken from the total coliform concentrations at station 
no. 31 (TC^]_) f result In calculated profiles within the bay which 
exceed observed levels. Knowing that station no, 31 is located in 
the ship channel (see Fig. 4.1) and that the concentration measured 



Table 4.3 Loadings of Total Conform at 


January 


February 


March 


April 


May 


Mouth of 
Mobile 
River 

Mouth of 
Tensaw 
River 

20,500 

2,000 

18,125 

2,000 

99,000 

2,000 

54,000 

2,000 

■ 

40,000 

200 

700 

300 

3,600 

1,000 


Mouth of 
Log 
River 


19,000 


13,800 


47,500 


7, 750 


1,800 


330 


330 


August 


1,500 


200 


200 


Alabama 

Port 


Cedar 

Point 


Mouth of 
Bon Secour 
River 







23,800 

2,500 

1,500 


5,000 

M50 

1,300 


2,100 

1,100 

170 


2,750 

550 

120 


1,100 

200 

40 


15 

1 

8 


60 

0 

45 






15 


2 


20 
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there is a point concentration instead of one characteristic of the 
entire model cell, a conversion of the point source concentration to 
a cell loading concentration suitable for model input becomes neces- 
sary. This conversion requires a dilution of the point source conc- 
entration to one distributed through the entire grid cell correspond- 
ing to station no. Jl. By definition, the dilution factor can be 
expressed as 


D.F. 


( volume of the actual water mass \ 

possessing the total coliform concentration ] 
observed in the field sample 


\ 


! 


volume of the grid cell corresponding to 
the location where field sample was taken 


However, due to the lack of detailed information about the magnitude 
of the waste discharge at point sources, and to the irregularity of 
the configuration of the water mass at the points of sampling, the 
dilution factor (D.F.) is determined by a calibration method involv- 
ing actual data. Point source data collected at coliform stations in 
the bay, especially samples collected in the ship channel or near 
possible waste outfalls where non-homogeneous mixing may exist, may 
not be representative of cell concentrations utilised in the NCSTM. 
Care must exercised in interpreting such kinds of data prior to their 
use as model input or for comparison purposes. 

Model Calibration 

The calibration procedure involved the adjustment of the 
source loading dilution factor (D.F.) and the adjustment of the 



7 ? 


correction factor (C.F.) for the dispersion coefficients according 
to Equation (3-40 ). A trial and error method was used, based on the 
June, 1962 total coliform data in which both the D.F. for point source 
loading concentrations and the C.F. for dispersion coefficients were 
varied. Values of these factors producing most consistent results 
over the entire range of the calibration data were selected and fixed 
for final use in the verification and the parametric phases of this 
project. The reason for choosing June, 19&2 data for model calibrat- 
ion was that for this month the river flow rate of the Mobile River- 
Tensaw River system (l) is close to the average value of 59i000cfs, an ^- 
(2) corresponds closely to the acceptable verification levels of the 
Hydrodynamic Model reported by Hill and April v The river flow 

rate of May, 1962 is closer to 59 ,000 cfs; however, the total coliform 
data for May is not satisfactory for verification purposes (see Table 

4.1 and Figures 4.11 to 4,23)* 

4.2 Results of the Verification Study 

The results are shown tabulated in Tables 4.4 to 4.11 for each 
month during which the verification phase of this study was conducted. 
Included in each table are the monthly mean total coliform concentrat- 
ion (in MPN/l00ml) and the 70 % confidence ranges calculated for the 
field data on a monthly basis. These values are compared with model- 
predicted total coliform concentrations for the cells corresponding 
to the stations in the bay where point concentration data are availa- 
ble. Furthermore, Figures 4.3 to 4.10 show the model- calculated total 
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coliform profiles within Motile Bay for each month from January to 
August , 1962, during which the verification phase is performed. 

Total coliform concentration vs. time (month) curves are also presen- 
ted to indicate the trend of concentration changes with season. Data 
at some stations lying very close to the tay perimeter were not 
selected for comparison, since data at these stations are not repre- 
sentative of the cell concentrations. 

Notations for Tables 4.4 to 4.11 

x - monthly average total coliform field data, 
t = the statistic used for confidence range correlation 
= S/frT = stadard deviation of the mean 

where S - standard deviation of the data 

n - no. of field samples for the month 
TC^p = total coliform concentration field data at station no. Jl 
(corresponding to the mouth of Mobile River) 

D.F. = dilution factor to convert TC^p into a cell loading source 
concentration 

=dieoff rate constant of total coliform bacteria 
E = dispersion coefficient 

Notations for Figures 4.11 to 4,23 

Actual total coliform concentration data 

Mo del -calculated monthly average total 
coliform concentrations 

Less than 
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Table 4.4 Total Coliform Concentration for Mobile Bay - January 1962 
Loading at Mobile River Mouth = l/4 TG^i ( B.F. -■ 4 ) 
Correction Factor for 3 ™ 500 


K r ~ 0.26 d.ay~^ 




* • 

• ♦ 


.ated Total Coliform Concentration 
* January, i 962 
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Table 4,5 Total Coliform Concentration for Mobile Bay -Februaryl962 
Loading at Mobile River Mouth — 1/5 TCaj_ ( D.F. - 5 ) 

Correction Factor for E ~ 500 
Kp = 0.29 day - ^ 


r — 

I 

! 

Station 
No , 

Measured Data 

1 

! 

j Model 

i ! 

I Calculated i 

1 

Result j 

No, of 
Field 
Sampling 

1 

Monthly 

Average 

X 

— 

70 % Confidence Range 

3 ? - tSj 

x + tS^ 

33 

2 

4,500 

CO 

<r\ 

r\ 

8,662 

1,981 j 

L 36 . 

4 

23,000 

18,125 

27,875 

< 

15,415 

1 

r 

| 39 

4 

5,000 

3,312 

6, 688 

i 

5,273 

! 60 

4 

17,000 

10,750 

23,250 

i 

9,677 : 


4 

63.500 

41 ,000 

86,000 

! 

10,548 

j 62 

j °5 

4 

27,500 

13,750 

41,250 

10,549 

4 

1,650 

881 

2,419 

3,514 

66 

4 

8,000 

7,531 

8,468 

6,421 

67 

4 

51,500 

33,562 

69,438 

8,407 

75 

4 

15,000 

7,312 

22,688 

j 

5,302 i 

i 

83 

3 

1,100 

300 

f 

1,900 | 

i 

^ 

1,491 

88 

4 

5,300 

1,800 

8,800 

3,783 

112 



3 

1,380 j 

603 

2,156 

709 
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Table 4.6 Total Coliform Concentration for Mobile Bay -March 1962 
Loading at Mobile River Mouth = i/5 TC^ ( D.F. - 5 ) 
Correction Factor for E = 500 
% = 0.39 day' 1 


Station 

No. 

Measured Data 

. 

| 

Model 

Calculated 

Result 

No. of 
Field 
Sampling 

Monthly 

Average 

* 

70$ Confidence Range 

X - tSjf 

x + ts* 

1 

O'N ! \0 

C°4 | Cn 

2 

8,000 

4,126 

11,874 

3,938 

3 

25,000 

13,638 

36,362 

80,734 

59 

4 

160,000 

-27,500 

347,500 

20,108 

60 

4 

69,500 

49, 688 

89,313 

41,863 

6l 

4 

35,000 

15,000 

55,000 

43,519 

62 

4 

14,000 

6 , 625 

21,375 

35,070 

65 

4 

4,160 

4,060 

4,260 

11,815 



66 

4 

36,000 

11,625 

60,375 

25,338 

67 

4 

19,250 

8,625 

29,875 

31,766 

75 

4 

15,750 

1,763 

29, 738 

17,326 

83 

4 

255 

186 

324 

3,159 

68 

4 

2,800 

1,363 

4,283 

5,3?5 

112 

3 

55 

-3 

113 

1,089 



























































Table 4.7 Total Coliform Concentration for Mobile Bay - April 1962 
Loading at Mobile River Mouth ~ l/4 TC^ ( L.F. ~ 4 ) 


Correction Factor for E = 500 
K r = 0.50 day“* 




Measured 

Data 



r | 

Model 

Station 





No . 

No. of 

Monthly 

70% Confidence Range j 

Calculated 


Field 

Average 



Result 


Sampling 

X 

1? - tSx 

x + tSjg 



! 

33 

: 

5 

. - — — 

1,540 

1,008 

2,072 

- — i 

2,727 ! 

i 

36 

5 

76,600 

44 , 669 

108,531 

44,613 
■ - — 

59 

4 

162,000 

-50,500 

374,500 

12,425 

60 

— 

4 

7,250 

5,688 

8,813 

24,451 

61 J 

4 

27,500 

14,313 

40,688 

23,589 

62 

1 

4 

17,000 

7,000 

27,000 

1 

9,722 

65 

1 

4 

7,100 

4,263 

9,938 

7,201 

66 

4 

8,100 

2,725 

13,475 

15,166 

67 

4 

2,750 

2,063 

3,438 

16,592 

75 

4 

5,600 

1,975 

9,225 

9,283 

! 83 

4 

30 

25 

35 

2,094 

j 

1 88 

. 

5 

850 

488 j 

1,212 

3,349 

112 

4 

55 

44 j 

66 

638 


























































o 


,ed Total Coliform Concentration 
.pril, 1962 
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Table 4,8 Total Coliform Concentration for Mobile Bay - May 1962 
Loading at Mobile River Mouth = 1/5 TC^} ( D.F. * 5 ) 
Correction Factor for E = 500 
K r = 0.72 day" 1 


Station 

No. 


33 

36 


59 


60 


61 


62 


66 


67 


75 


83 


88 


112 


No. of 
Field 
Sampling 


Measured Data 


Monthly 

Average 

x 


150 


91,600 


600 


10,600 


6,000 


670 


3,500 


5,240 


20,000 


3,000 


19 


260 


12 


70% Confidence Range 


X - ts x 


56 


25,07? 


-13 


-1,108 


-1,108 


91 


906 


717 


5,578 


925 


13 


x + tS x 


244 


158,123 


16 


1,213 


22,308 


13,108 


1,249 


6,094 


9,763 


34,422 


5,075 


25 


504 


17 


Model 

Calculated 

Result 


1,515 


1,456 


3,287 


2,824 


1,250 


638 


1,523 


1,514 


582 


72 


434 


22 


1 




























































Fig. 4.7 Model- Calculated Total Coliform Concentration 
Profiles for May, 1962 




Table 4,9 Total Coliform Concentration for Mobile Bay - June 1962 
Loading at Mobile River Mouth = i/4 TC^ ( D.F. - 4 ) 
Correction Factor for E = 500 


K,= 0. 


Station 

No. 

Measured Data 

Model 

Calculated 

Result 

No. of 
Field 
Sampling 

Monthly 

Average 

X 

70% Confidence Range 

Y - tSx 

x + tSx 

33 

3 

250 

145 

355 

70 

36 

3 

600 

320 

880 

293 

59 

4 

12 

7 

1 ? 

26 

60 

4 

25 

10 

40 

61 

61 

4 

50 

27 

73 

81 

62 

4 

110 

19 

201 

124 

\0 

4 

20 

7 

33 

12 

66 

3 

7 

4 

10 

26 

67 

3 

132 

. ... . 

14 

250 

66 

75 

3 

20 

8 

32 

16 



83 

4 

10 

9 

11 

H 

88 

3 

15 

14 



16 

7 

112 

4 

10 

6 

14 

1 

















































































Table 4.10 Total Coliforjn Concentration for Mobile Bay - July 19°2 
Loading at Mobile River Mouth - l/6 TC-^ ( D.F. - 6) 

Correction Factor for E = 500 
K r = 0.88 day" 1 


Station 

No. 

Measured Data 

Model 

Calculated 

Result 

No. of 
Field 
Sampling 

Monthly 

Average 

X 

?0$ Confidence Range 

X - tSx 

x + tSjj 

33 

4 

360 

154 

566 

249 

36 

4 

300 

144 

456 

1,272 



59 


9 

6 

12 

86 

60 

i 

4 

35 

12 

58 

176 

6l 

4 

161 

69 

101 

166 

62 

4 

100 

50 

150 

138 

65 

4 

20 

9 

31 

40 

66 

4 

40 

15 

65 

78 

6? 

4 

33 

13 

53 

86 

75 

4 

120 

10 

230 

29 

83 

5 

10 

9 

11 

3 

88 

5 

13 

5 

21 

20 

112 

5 

40 

20 

60 

2 
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Table 4.11 Total Coliform Concentration for Mobile Bay -August 1962 
Loading at Mobile River Mouth - i/6 TC^ ( D.F. ~ 6 ) 

Correction Factor for E = 500 
K r = 0.90 day - * 


Station 

No. 

Measured Data 

Model 

Calculated 

Result 

No. of 
Field 
Sampling 

Monthly 

Average 

X 

70% Confidence Range 

X - ts x 

x + tS x 

33 

5 

50 

32 

68 

74 

36 

5 

250 

160 

340 

528 

59 

2 

4 

2 

6 

36 

60 

2 

10 

8 

12 

104 

61 

2 

15 

1 

31 

162 

62 

3 

15 

12 

18 

282 

65 

3 

10 

3 

17 

15 

66 

3 

7 

4 

10 

40 

6? 

3 

5 

4 

6 

89 

75 

4 

4 

3 

5 

27 

. — 

83 

4 

8 

5 

11 

1 

88 

3 

3 

2 

4 

8 

112 

4 

7 

1 

13 

1 
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Fig. ^.17 Model calculated averages compared with actual data of total coliform concentration 
at station No. 65 « 
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Fig* 4.20 Model calculated averages compared with actual data of total coliform concentration 
at station No. 75* 
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Fig. 4.22 Model calculated averages compared with actual data of total coliform concentration 
at station No. 88. 
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4. 3 Discussion of Verification Kesults 

Verification of the NCSTM for total coliform is based on the 
comparison between the model -calculated monthly average results and 
the monthly average actual field data. The total coliform concentr- 
ation, unlike DO levels which usually become saturated between 7 to 
15 ppm (see Fig, 3*4 on p. 4?) and estuarine BOD levels, which 
normaly are below 1000 ppm, is a much more variable entity described 
by rather broad ranges of numerical values. This results in the 
restriction of the verification phase to essentially trend-analysis 
levels. 

Numbers of field samples varied from 2 to 5 per month for each 
station during the months January to August, 1962 (see Table 4.4 to 
Table 4.1l). These data are scattered as a result of varying field 
conditions and sampling accuracy, and when averaged over the monthly 
periods, wide variations in the standard deviations occur. The 
standard deviations, of these samples were used to calculate the 
standard deviations of the mean, S^, for each datum on a monthly 
basis. Sx was in turn used with the t-distribution to calculate the 
70$ confidence range for each datum for each month. Details of the 
use of the statistic "t” are cited by Volk^O). 

Trend Analysis Verification Comparisons 

The verification comparison (trend analysis) consists of three 
steps, as summarized below. 
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Step 1* Ihe comparisons between the model -calculated results and 
actual field data for each of the 13 stations shown in 
Figures 4*11 to 4.23 indicate that the model is capable of 
following the trends of the total coliform concentration 
within Mobile Bay. Except for station no. 59 (see Fig. 4.5)» 
where extraordinarily high monthly mean concentrations were 
measured for March and April, most of the trends for season- 
al variation shown for other stations are reasonably accura- 
te. Stations 36 , 60, and 61 are located in the chip channel 
(see Fig. 4.l). Deviations in model -predicted results for 
these stations are found to be in agreement with expected 
trends based on the hydrodynamic behavior of the ship 
channel. Station no. 88 is located near Cedar Point; 
stations no. 83 and 112 are located in Bon Secour Bay (see 
Fig. 4.l). These three stations tend to show deviations 
which are more pronounced in these regions as a result of 
seawater intrusion from the Gulf of Mexico. It is found that 
except for station no. 112 in February, all the deviations 
can be explained by the seawater intrusion process; that is, 
positive deviations (calculated results greater than actual 
data) in cold months during which actual K r values are lower 
in the bay area than in the Gulf, and negative deviations 
(calculated results less than actual data) in warm months 
during which K r values are higher in the bay area than in the 
Gulf, The reflectional month is May at which time total 
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coliform concentrations at all stations change drastically. 
Total coliform concentration data for many stations for the 
periods October to December, 1962 and September, 196I are 
not available for the water year 1962, and therefore they 
are not included in this study. 

Step 2i By inspection of the total coliform concentration profiles 

shown in Figures 4.3 to 4.10, the western half of Mobile Bay 
is usually suffering more severely from pollution than the 
eastern half of the bay. At the same time during the months 
of high river flow rates, there are very obvious "tongue” 
effects in which the total coliform concentration profiles 
reach far down the bay. During months of low river flow 
rates this "tongue” effect becomes much smaller* These 
behaviors are coincident with what has been evaluated ir the 
work by Gallagher, et_al, and shown by the studies of 
Hill and April ( l2 ). 

Step 3t By inspection of Figures 4.11 to 4.23, it is found that 

model- calculated results fall within the ranges covered by 
the actual field data 66 , % of the time for the period 
January to August, 1962. It should be indicated that this 
percentage includes stations located within the ship channel 
which are not expected to show good agreement with model - 
calculated results as a consequence of the way in which 
the model is formulated (i.e. two-dimensional, no stratifi- 
cations). Due to the inclusion of the ship channel data in 
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addition to the obvious scatter in the available field data, 
model verification for total coliform can only be made for 
trend analysis purposes. More detailed point analyses must 
be deferred to a time when more accurate and precise field 
data measurements can be obtained. 

More specific factors relating to the verification of the 
model -predicted total coliform distribution are itemized and their 
effects are discussed in the following paragraphs. 

Dilution Factor Correlation 

Dilution factors (D.F.) were varied between 2 and 10 in the 
preliminary model calibration study. A D.F. of 4, together with a 
value of 500 for C.F. for the dispersion coefficients (see Eq_. 3“40 
on p. 39) were found to produce reasonable total coliform distribu- 
tions based on the June, 1962 data. This set of factors was fixed to 
exercise the NCSTM for other periods for verification purposes. 
Subsequent fine tuning of the model in the verification phase showed 
that the D.F. may be regarded as a function of the river discharge 
rates and the total coliform source concentration at the mouth of 
Mobile River draining into the bay. Different river discharge rates 
result in different flow velocities and different degrees of mixing. 
In a portion of the verification study refinements of the dilution 
factor values indicated the followings a dilution factor of 4 was 
required for the months with highest river discharge rates (January 
and April), a dilution factor of 5 was required for months with 
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medium river flow rates (February and March), and a dilution factor 
of 6 was required for months with low river flow rates (July and 
August) . This is in agreement with the fact that mixing is greater 
for high river discharge rates and therefore a smaller dilution factor 
would be required (i.e. the grab sample would be more closely repres- 
entative of the cell concentration.) Verification results for June 
are exceptions to this trend, where a dilution factor of 4 was used 
with good results. Note, however, that the source loading concentra- 
tion for this month at the mouth of the Mobile River (i.e., TC^l) was 
extremely low compared to other summer months (see Table 4,3)» while 
the river discharge rate remained high. This may be explained by the 
fact that better mixing was attained at station no. 31, and therefore 
a small dilution factor was required for the conversion to a cell 
input concentration, A D.F. of 4 used for May gave best results 
compared to the monthly averages of data collected for that month. 

The total coliform group concentration for all stations within the 
bay undergoes nearly a step change from high levels to low levels in 
May, Data are especially scattered, as shown in Table 4.8, and the 
comparison should be regarded as less significant. 

For points where serious pollutant transport is expected, 
such as at points near the mouths of the rivers and loading sites, 
calculated values of total coliform concentration are usually smaller 
than the actual data. This is because the model calculates cell 
concentrations, while actual data are grab samples collected from 
water which is not well mixed. Hie ship channel is rather narrow as 



113 


compared to the grid cells that encompass it. Concentrations calcu- 
lated for grids located in the ship channel are usually smaller than 
the actual data; that is, at stations no. 36 » 60, and 66, etc. This 
is due to the high current velocities in the ship channel; the reten- 
tion time for total coliform "bacteria is reduced, resulting in higher 
point concentrations compared with adjacent, slower moving bay waters. 
For higher river flow months this effect was so pronounced that total 
coliform profile contours reached far down the left half of the bay 
(see Figures 4.4 to 4.7). This is consistent with the observations 
by Gallagher, et al . and Hill and April^"^. 

Dispersion Coefficients 

Dispersion coefficients are calculated by Eq. (3-40). Values 
of 250, 500, 750 t and 1000 were tried for the correction factors (C.F.) 
on both x- and y- component maximum tidal velocities in the model 
calibration phase. It was found that a value of 500 gave reasonably 
good results for both the x-component and the y- component dispersion 
coefficients, Ex and Ey. It was also found that smaller C.F. values 
usually elevate coliform concentrations in the upper portion of the 
bay near the waste sources and decrease coliform concentrations in 
the portions of the bay far removed from .the waste sources . For each 
monthly period, the smaller the C.F. used, the more pronounced is this 


observation. 
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Dieoff Rate Constants 

Larger values of the dieoff rate constant K r decrease the 
total coliform concentration values in any given location id. thin the 
"bay. The opposite is true for smaller K r values. Therefore it is 
important that correct temperatures be used to calculate the corres- 
ponding K r values needed by the model equation. 

Seawater from the Gulf of Mexico and from the Mississippi 
Sound causes a slight temperature shift within the southern section 
of the bay. Seasonal average sea surface temperatures of the Gulf 
of Mexico, obtained from the National Atlas of the U« S. A.(^), are 
listed in Table 4.12. Temperatures used in the verification are also 
listed for comparison. It is seen that the temperature of Gulf waters 
is more stable, i.e., is varying over a smaller range than that 
of the bay water. Due to seawater intrusion, the water temperatures 
in the lower portion are also more stable than those at other portions 
of the bay. Since dieoff rate constants are directly related to 
temperature, it is expected that actual concentrations (data) at those 
stations in the lower portion of the bay will be affected. This is 
particularly true in the Bon Secour Bay area and for points near 
Main Pass and Cedar Point. Observed concentrations should be lower 
in warm seasons and higher in cool seasons, as compared to what would 
be calculated by the model based on homogeneous bay temperatures. 
Station no. 88 is located near Cedar Point j station no. 83 and station 
no. 112 are located in the Bon Secour Bay area (see Fig, 4.1). They 
are subjected to intrusion of seawater from the Mississippi Sound and 
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Tb/ble 4.12 Comparison of the surface water temperature of the 
Gulf of Mexico near Motile Bay and the bay water 
temperatures used in verification of the NCSTM for 
total coliform. 


Seasons 

Gulf surface temperatures, 

Op 

Bay water temperature 

Maximum 

Minimum 

Average 

used for verification 
study, °F 

Spring 

82 

70 

76 

March 6l . 3 
April 67*9 
May 78.1 

Summer 

. 

86 

78 

82 

June 81.4 
July 83.7 
August 84. 2 

Fall 

78 

64 

71 

September 

October 

November 

Winter 

74 

58 

66 

December 

January 49 . 5 

February 58.3 




























116 


the Gulf of Mexico. Most data at these stations are noted to be 
different from calculated results in the described direction. The 
extent of seawater intrusion is seasonal in nature, and therefore the 
area affected by this phenomenon varies. In this study it is found 
that except for stations no. 83 and 112 in the months of March and 
April, all the deviations are small (see Tables 4.4 to 4.1l), and for 
the purpose of trend analysis based on monthly average calculation, 
this effect can be neglected without introducing too significant 
errors . 

In the following chapter, attention will be turned to the 
parametric study in which the sensitivity of the mo del -predicted 
results toward various changes in system behavior was investigated. 



CHAPTER V 


PARAMETRIC STUDY 

There are four major parameters which affect the total 
coliform distribution in Mobile Bay. These are* 

(1) river flow rates, which influence the total coliform concentration 
introduced into the bay and the retention time of the bacteria within 
the bay; 

(2) wind conditions, which influence the current distribution and 
therefore the retention time of bacteria within certain portions of 
the bay; 

(3) temperature, which influences the death rate of total coliform 
bacteria; and 

(4) waste loadings, which influence the input concentration of total 
coliform bacteria introduced into the bay from various sources. 

These variables are examined in a parametric study to determine how 
sensitive the total coliform group concentrations are to changes in 
variable magnitude and/or direction which simulate real system condit- 
ions. The first three of the four parameters, i.e., river discharge 
rates, wind conditions and temperature, will be studied in Section 5«1» 
the waste loading effect will be studied in Section 5*2. 
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5.1 Effects of River Flow Bates, Wind Conditions, and Temperature 

Table 5*1 shows the input data used for the 18 parametric 
runs performed to study the effects of river flow rates, wind condit- 
ions, and temperatures. Two levels of wind, i.e., 15 knots and 25 
knots, are studied and compared with results calculated for no wind. 
Three directions of wind are studied. They are: from the north 
( 0 = 90° ), from the southwest ( 0 = 225° ), and from the southeast 
( 0 - 315° )• The value 0 is the wind direction in degrees, as mea- 
sured in the counter-clockwise direction from the x-axis in the 
Cartesian -coordinate system. For medium river flow rates all three 
directions are studied to determine the effect of variation of wind 
direction (see Table 5*1 > runs a to g). For low and high river flow 
rates (see Table 5»1» runs h to j and k to m), wind from the southwest, 
the most prevailing direction, is studied. 

Total coliform source concentrations used in the parametric 
study come from the data used in the verification analysis (1962 
period) having comparative levels of river flow rates. Thus those 
loading concentrations of May, 1962 are used for medium river flow, 
those of August, 1962 are used for low river flow, and those of April, 
1962 are used for high river flow (see Table 4,3 on p, 75), to exerci- 
se the model. 

To study the effect of variations of river flow rates, the 
conditions experienced in May, 1962 are used as a reference. The 
reason for selecting this month is that river flow rates are more 
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Table 5*1 Data Used for Parametric Study Puns a to r 


Run 

Wind 

River Discharge Rates 

Temperature 

Speed 

knots 

e 

deg. 

Mobile 

River 

cfs 

Dog 

River 

cfs 

Tens aw 
River 

cfs 

°F 

a 

0 

- 

24,000 

2,000 

20,000 

78.1 

b 

15 

225 

II 

11 

tl 

M 

c 

25 

225 

It 

11 

fl 

« 

d 

15 

90 

ff 

ft 

ft 

if 

e 

25 

90 

11 

ft 

• t 

tt 

f 

15 

315 

it 

M 

II 

tt 

6 

25 

315 

ft 

. 

ff 

ft 

h 

0 

- 

7, 000 

500 

5,000 

84.2 

i 

15 

225 

If 

If 

i» 

IP 

j 

25 

225 

11 

if 

n 

ft 

H 

0 


145,000 

5,000 

100,000 

67.9 


15 

225 

ft 

II 

it 

it 

m 

25 

225 

fl 

ff 

ft 

ff 

n 

7.9 

225 

10,000 

1,000 

9,250 

78.1 

0 

It 

n 

40,000 

4,000 

37,000 

9f 

P 

fl 

tv 

20,000 

2,000 

18,500 

PI 


If 

«* 

ff 

ff 

ti 

85.8 

r 

U 

it 

ft 

ff 

tt 

68.1 


Die- 

off 

Rate 

K 

day - ' 

0.72 

»l 

11 


If 

ff 

0.90 

ft 

If 

O.50 

If 

ft 

0.72 

II 

I* 

0.94 

0.50 
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representative of average conditions and those conditions under which 
the Hydrodynamic Model for Motile Bay was initially verified (see 
Table 4.1 on p. ?0). By holding wind conditions and temperature 
constant, the value of river flow rates is first doubled, then halved 
(see Table runs n, o, and p), to determine the effect on total 

coliform concentration within the bay. 

Similarly, to study the effect of variations of temperature, 
May, 1962 data, are again used as a baseline in which river flow rates 
and wind conditions are held constant. The temperature used to exer- 
cise the model is first raised from ?8,1°F to 85.8°F (see Table 5«1» 
run q) to give an increase in the dieoff rate constant K from 0.72 
day" to 0.84 day . The temperature is then reduced to 68.1°F (see 
Table 5-l> 3?un r) to give a decrease in the dieoff rate constant from 
0.72 day - ^ to 0.50 day”^ to exercise the model. 

The computational procedures are similar to that used in 
the verification study (see Sections 3 . 4.3 and 4.1), The results of 
the parametric study runs as listed in Table 5.1 are shown in Figures 
5.1 to 5*9* On each figure, comparisons are made at three levels of 
total coliform concentration in units of MPN/lOOml. Each coliform 
concentration contour is marked with the letter identifying the 
corresponding parametric study run. The way parametric runs listed 
in Table 5»1 are combined for various comparison purposes is given in 
Table 5-2. In the following paragraphs, discussions of the effect 
each variable has on total coliform distribution within Mobile Bay 
are presented. 
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Table 5.2 List of Figures for Parametric Study Comparisons 


Figure 

comparison 

among 

runs 

indicating 
the effect of 
variation of 

at 

constant 

5-1 

'p q. r 

temperature 

river flow rates and 
wind 

5.2 

nop 

. ... 

river flow rate 

wind and temperature 

5.3 

a b c 

speed of wind from SW 

medium river flow 

5.4 

a d e 

speed of wind from N 

medium river flow 

5.5 

a f g 

speed of wind from SE 

medium river flow 

5-6 

h i j 

speed of wind from SW 

low river flow 

5*7 

k 1 m 

speed of wind from SW 

high river flow 

5.8 

a b d f 

direction of wind at 
15 knots 

medium river flow 






5-9 


a c e g 


direction of wind at 
25 knots 


medium river flow 






































,tion profiles from runs 
.ts of profiles due to 

r rates 









.ration profiles from runs 
river flow rates j wind from 

















,ion profiles from runs 
l river flow rates? wind 
■om SW, N, and SE. 
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Temperature 

Fig. 5.1 shows the effects of changing temperatures on total 
coliform distribution. The shifts of the 100 and the 5^0 MPN/lOOml 
total coliform concentration isolines as seen in Fig. 5»1 are 5- n the 
order of 2 to 4 grid widths (4 to 8 km.) from run to run, which can 
seriously affect the shellfish harvesting activities in the bay, 
especially in the Bon Secour Bay area. This simulates what can 
happen to the coliform distribution in case of sharp temperature 
variation, when all the other system variables, i.e., river flow rates, 
wind conditions, and waste loadings, remain unchanged. The reason for 
such pronounced shifts of coliform concentration profiles is the 
change in dieoff rate constant, K, caused by temperature variation. 

The change in K follows Eq. ( 3 -^ 2), which indicates that K is a funct- 
ion of bay water temperature alone. When water temperature in the bay 
is higher, total coliform bacteria dissipate at a higher rate, and the 
coliform concentration in the bay becomes lower. When the water 
temperature is lower, K is smaller, the total coliform bacteria die 
off at a lower rate, and the coliform concentration in the bay becomes 
higher. This effect also partly accounts for seasonal variation of 
total coliform concentration within Mobile Bay. Seasonal variations 
of coliform concentration profiles have also been depicted in the work 
of Gallagher, et al . 
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River Discharge Rates 

The effect of variations in river discharge rates on the 
total coliform distribution profiles in Mobile Bay is shown in 
Fig. 2. The values of the river discharge rates using May, 1962 
data as a baseline are first doubled, than halved, to run the model. 
The results are then compared to the actual flow rate observed during 
1962. For decreases in river flow rates, the contours obviously shift 
upward, which results in a lower overall coliform distribution within 
the bay. For higher river flow rates, the contours all shift down- 
ward, which results in a higher overall coliform distribution. The 
reason for these changes is two -fold. By holding the loading conc- 
entrations constant and increasing the river flow rates, more total 
coliform bacteria are introduced into the bay, while at the same time 
the net current velocities in the negative y -direction (north to south) 
are increased. This latter condition allows less retention time for 
the total coliform group to die off, and results in higher residual 
coliform concentrations at any part within the bay. For lower river 
flow rates the reverse is true. These effects caused the changes 
observed in the runs in Fig. 5 . 8 , and are consistent with actual 
observations in Mobile Bay^ 0 ^. 

When river flow rates are higher due to either rainfall or 
storm, the amount of coliform group bacteria loaded into the river 
water by runoff is indeed higher. However, the loading concentration 
at those loading grids may not be constant. In parametric runs n, o, 
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and p they axe assumed to "be constant. This assumes a linear relation 
"between amount of total coliform group input and fresh water runoff, 
and neglects the difference in fresh water runoff from agricultural 
areas and those from municipal areas. A more realistic way of assess- 
ing the effects of changing waste loading independent of river dis- 
charge rate is discussed in Section 5-2. 

Wind Effect 

For medium river flow rates (see Table 5*l)» Fig. 5*3 to 5-5 
show the effects of variations in wind speed (0, 15 » and 25 knots) 
blowing from three different directions (N, SE, SW). Fig. 5*6 and 
5.? illustrate for low and high river flow rates respectively, the 
effects of changing wind speed originating from the southwest directi- 
on. Fig. 5*8 and 5-9 show the effects of variation in direction of 
wind at 15 and 25 knots, respectively, at medium river flow rates. 

In each comparison, the temperature of the bay water is held 
constant, and the observed variation is exclusively due to variation 
in net velocities and dispersion coefficients resulting from varying 
wind and river discharge values. Increasing the net current veloci- 
ties in the negative y-direction (from the north to the south) will 
shorten the retention time the coliform bacteria would spend within 
the bay, allow less time for coliform to die off, and thus increase 
the total coliform concentration at any location within the bay. From 
the model calibration study of Chapter III it has been found that 
deliberately increasing the dispersion coefficients (by manipulating 



the correction factor for dispersion coefficients) would tend to 
decrease higher coliform concentrations and increase lower coliform 
concentrations, i.e., would tend to shorten the range of concentrati- 
ons. In the parametric study comparisons (Fig. 5*3 to Fig. 5 * 9 ) » 
however, the change in either net velocities or dispersion coefficie- 
nts is neither uniformly increasing nor uniformly decreasing for all 
the grid cells. The final changes in coliform distribution are the 
gross totals of the effects of all the local changes in x- and y- 
component net velocities and dispersion coefficients. For most 
comparisons it may "be summarized that the displacements of profile 
contours are in the direction of the wind, i.e., the winds have 
caused the profiles to shift in the directions of the winds. However, 
the displacements rarely exceed the width of one grid (2 km). Run 
identifications are used to indicate their relative positions, as 
shown in Figures 5*3 to 5 * 9 * For many cases, the profile contours 
are so close together that their difference are not discernible. Due 
to the fluctuations in physical environments, these displacements 
would readily be masked and become undetectable. Therefore, the eff- 
ects of changing wind speed and direction on the monthly average 
total coliform distribution within Mobile Bay can for all practical 
purposes be regarded as negligible. To determine if wind has a great- 
er influence on distribution of total coliform for periods less than 
one month, more detailed data must be used to verify the model. These 
data are not available at the present time. 
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5.2 Effect of Waste Loadings 

Cell loading concentration of total coliform at the mouth of 
a river reflects the pathogenic pollution potential the river has on 
the bay. This concentration is contributed by waste loadings from 
various sources such as municipal, industrial, and rural areas. 

Table 5»3 shows the parametric runs performed in the study on the 
effect of changing levels of waste loading, which are expressed in 
the form of total coliform concentrations at the boundary cells 
representing the mouths of the rivers. The conditions experienced in 
May, 1962 are again used as a reference. River flow rates, wind 
conditions, and temperature are held constant. The only changes made 
are on the loading concentrations of total coliform bacteria at the 
mouths of Mobile River and Dog River. Values are reduced to l/2, l/4, 
and l/8 of the values experienced in May, 1962 to exercise the NCSTM. 
The resulting total coliform concentration profiles are shown in Fig. 
5.10. Comparisons are made at two concentration levels, i.e,, 70 and 
1000 MPN/lOOml. Ekch concentration contour is labeled with the letter 
identifying the corresponding parametric study run. Fig. 5*10 shows 
that each of the shifts of the coliform concentration profile is in 
the order of 2 grid widths (4 km). It is noted that the 70 MPN/l00ml 
contour shifts as many as 6 grid widths from run p to run u, as 7/8 
of the original total coliform bacteria is removed or reduced. These 
changes in total coliform loading are more realistc of conditions 
that might be achievable for varying degrees of treatment. 
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Table 5.3 Tata used for parametric runs p, s, t, u. 


Run 

River Flow Rates, 

; 

loading 

Concentration (MPN/lOOml) 
at 


Wind Conditions , 
and Temperature 

Mobile 

River 

Mouth 

Dog 

River 

Mouth 

Other 
lo cation 


Same as run p 



Same as those 

, 

in Table 5-1 

*+0,000 

1,800 

of May, 1962 in 
table 
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This study provides a method for describing total colifomt 
bacteria concentration distributions in Mobile Bay, and for describ- 
ing how these distributions might be affected by various changes in 
the real system. As a result of these preliminary investigations, 
a series of conclusions and recommendations related to the use and 
extension of ideas generated within this study is presented in 
Chapter VI. 



CHAPTER VI 


' CONCLUSIONS AND RECOMMENDATIONS 

Chapters IV and V presented results demonstrating the feasib- 
ility of using a two dimensional model to describe the transport of a 
non-conservative species within Mobile Bay. The particular species 
investigated was the total coliform bacteria group. The intent of 
Chapter VI is to present the concluding observations from this study, 
the limitations of the present model, the contributions resulting from 
this study, and the recommendations for continued research in related 
areas . 

Concluding Observations 

A model for the prediction of trend behavior of the total 
■coliform bacteria distribution within Mobile Bay has been developed. 
This model allows fors 

(1) variability in the total coliform source concentration at several 
locations along the boundary of the bay system, 

(2) variability of the correlation coefficients for the x- and 
y-component dispersion coefficients to best describe the mixing 
characteristic of the specific non-conservative species, 

( 3 ) variability in river flow and wind conditions by interacting with 
the Hydrodynamic Model developed by Hill and April 
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(4) variability in temperatures that result in changes in the dieoff 
rate constant of the total coliform bacteria. 

This model is based on established engineering practice and 
constitutes the necessary framework for the development of other 
similar non-conservative species transport models for BOD and DO. 
Additionally, the model formulation is made in such a manner as to 
facilitate rapid execution and easy interpretation of computed results, 
which are printed in the same configuration as the bay. 

Specific observations related to various phases of this study 
should also be presented. In the verification phase of this study, 
the model was calibrated with the June, 1962 total coliform concentr- 
ation data taken from the bay. Values of the dilution factors (for 
the conversions of point source loading concentrations to cell loading 
concentrations suitable for model input) and the correction factors 
(for the correlation of x- and y-component dispersion coefficients 
based on the x- and y-component maximum current velocities over the 
tidal cycle) that best describe total coliform mixing characteristics 
were calibrated. They were then used for the verification of the 
model based on the actual data collected during January to August, 

1962. In the parametric study phase of this investigation, tempera- 
ture was found to have the most pronounced effect on the total coliform 
distribution within Mobile Bay. A change in temperature of 10°F can 
cause the total coliform concentration profiles in the bay to be dis- 
placed as much as 8 kilometers. Variations in river flow rates also 
showed a pronounced effect upon the total coliform distribution pro 
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lies within the bay. The May, 19&2 condition was used as a baseline 
to study the effect of changing river flow rates. When the river 
flow rates were doubled and then halved to exercise the model, dis- 
placements of total coliform concentration profiles by as much as 
6 kilometers were obtained. Wind conditions (speed and direction) 
were studied at three speeds and three directions, interfacing with 
three levels of river flow rates. It has been found that wind condit- 
ions have the least influence on monthly average total coliform distr- 
ibutions within Mobile Bay as compared with other parameters. However, 
it is believed that reduction of the time basis to a tidal level will 
result in the observation of more pronounced wind effects than those 
observed from the monthly averaged results. 

In addition to the above, total coliform source concentration 
levels were varied at constant temperature, river flow rates, and wind 
conditions to simulate the possible effects different treatments would 
have on total coliform distribution within the bay. Conditions exper- 
ienced in May, 1962 were again used as a baseline. Source loading 
concentrations of total coliform experienced in May, 1962 were reduced 
to l/2, l/4, and l/8 to exercise the model. Displacements were found 
to be in the order of 2 to 8 kilometers from inn to run. 

Limitations of the Model 

At the present time, the most limiting factor involved with 
any modeling activity on Mobile Bay is the availability of suitable 
field data for calibration and verification of the formulated models. 
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This fact restricts the NCSTM developed in this study to a trend 
analysis tool. By this, it is meant that specific predictive capa- 
bilities related to source or cell concentrations can not be assured 
with any degree of confidence that the results reproduce real system 
behavior. However, monthly trend analyses of species concentrations 
for given regions of the bay system resulting from natural or man- 
made phenomena can be assessed with relatively high accuracy. 

Also, the present NCSTM is limited to those conditions for 
which calibration and verification of the Hydrodynamic Model were 
achieved. These conditions include (l) combined river flow rates 
of the Mobile River-Tens aw River system between 12,000 and 245,000 
cfs, and (2) wind speeds lower than 25 knots. Any conditions which 
approach the limits of the above should not be expected to produce 
reliable results unless further testing is made. 

Also included as limitations to the model formulated in 
this study are t 

(1) constant density of water throughout the bay, 

(2) normal tidal conditions at the Gulf boundaries, 

( 3 ) binary mixing behavior within each cell of the bay model, 

(4) homogeneous temperature throughout the bay for each 
month, and 

( 5 ) tidal average velocities and dispersion coefficients. 

These limitations should be reevaluated as more sophisti- 
cated model capabilities are developed and more reliable data are 
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obtained. 

Contributions of This Study 

The greatest contributions of this research lie in the 
development of a tool for the rapid assessment of conditions within 
Mobile Bay, and, the provision of a base from which other pertinent 
models may be developed. This study represents a continuing effort 
in the development of a comprehensive model for a detailed analysis 
of many proposed activities pertinent to a progressive society. 
Contributions were made, also, in areas summarized below. 

(1) A trend analysis of the total coliform was made, to better under- 
stand how this species is transported through the bay. This 
hopefully will lead to a better understanding of those variables 
affecting total coliform distributions, as that progress can be 
made to reduce their levels to allow for better economic growth 
within the shellfish harvesting industry. 

(2) The trend analysis should also provide insight into the deve- 
lopment of models for related non -conservative species which 

are indicators of water quality within the bay. These additional 
species include BOD and DO, which are widely accepted as standards 
for measuring industrial and municipal pollution loadings in 
natural water systems. 

( 3 ) The interactive effect of physical and biological terms has been 
demonstrated in this study, opening the door for interdiscipli- 
nary research and development projects. Through these inter- 
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disciplinary programs, better understanding of coastal ecosystems 
can be achieved. At a time when the coastal zones are being 
developed to provide the resources for energy related projects 
(i.e. deep water port development, off-shore and near-shore oil 
explorations, on-shore refining and processing facilities, etc), 
an understanding of the effects that these developments might 
have on this complex, interactive system is essential. 

Recommendations for Further Study 

Based on the experience gained in this study, several re- 
commendations are made concerning further studies in related areas. 

These recommendations are summarized below. 

(1) Establish a system within the bay area for routine, synoptic 
data collection in support of the mathematical modeling efforts. 
This could be achieved with little additional expenditures, pro- 
vided that cooperation among those agencies and organizations 
conducting active research programs within the bay can be esta- 
blished. 

(2) Investigate the use of the NCSTM to predict BOD and DO within 
the bay. 

( 3 ) Investigate more closely the mechanisms that govern the repro- 
duction, dieoff, and reactions of the various species related 
to water quality in the bay. 

(4) Identify the agencies within the region which have a need for the 
predictive capabilities of such models. These agencies include, 
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but are not limited to* 

Alabama State Department of Health 
Alabama Water Improvement Commission 
Alabama. State Geological Survey 
U. S. Corps of Engineers (Mobile District) 

Marine Environmental Sciences Consortium 
Alabama State Department of Conservation 
Private Industries including: 

Petroleum Companies, 

Petrochemical Companies, 

Chemical Companies, etc, 

which have active interests in developments within and adjacent 
to the bay. 

Hopefully these recommendations can be implemented to such 
an extent that mathematical modeling efforts can be used to assist 
in the development of protective systems for our coastal environment 
to keep pace with the ever-increasing population and industrial 
development. 
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A1 PROCEDURES FOR EXERCISING THE NCSTM 


(i) Elow Diagram Representation of the Procedures of 
Exercising the NCSTM for Mobile Bay 
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(il) Explanation 

River flow rates, wind conditions, hydrographic information, 
and other data were input to the Hydrodynamic Model for Mobile Bay 
developed by Hill and April For details of exercising the 
Hydrodynamic Model which is the first step of exercising the NCSTM, 
users are referred to the work by Hill and April. For the computation 
involving different periods, the main changes in input data are the 
river flow rates and the wind conditions. These changes can be made 
by altering two data cards. Two consecutive runs (the second one 
is initiated with the results obtained by the first one) are usually 
required to obtain convergent results of net velocities and maximum 
current velocities over the tidal cycle. These results are then 
stored into data files in the mass storage of the computer system, 
later to be retrieved by the NCSTM as input data. 

Bay water temperature, total coliform cell loading concentra- 
tions, and other pertinent data are input to the NCSTM together with 
the data created by the Hydrodynamic Model. The NCSTM then computes 
the desired total coliform distribution within Mobile Bay. The dis- 
tributions are printed with variable format to simulate the con- 
figuration of the bay . .Only total coliform concentration values at 
water cells are printed. Values at those cells representing the river 
channels and at those cells corresponding to land cells are not 
printed (values at land cells are zeroes). 

Computations were performed with a UNIVAC-1110 computer system 
on which provisions exist for users to store the main program in the 
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mass storage of the system. This saves the users the time required 
to read the source program deck into the system, especially when a 
large number of punched cards are used. In this study both the main 
program of the Hydrodynamic Model and that of the NCSTM involve large 
numbers of cards, and therefore are stored in the system. Descrip- 
tions of the control cards for the storage and retrieval of stored 
programs and calculated data are cited in Appendix A3. 
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A2 DESCRIPTION OP VARIABLES APPEARING IN THE COMPUTER PROGRAM 


Variable Name 
A 
At 

A2 

A3 

A4 

A5 

A6 

A? 

BRK 

C 

COLC 
. COLG 
COLM 
COLT 
D 

DELS 

BELT 

EX 

EY 


Definition 

= 6561 .68 ft 

storage term in the computation of cell total 
coliform concentration, in (MPN/lOO ml)* (ft/sec) 


storage term in the computation of cell total 
coliform concentration, in (MPN/lOO ml) 

dieoff rate constant of total coliform bacteria, 
(day" 1 ) 

total coliform concentration (MPN/lOO ml) 

total coliform cell loading concentration at 
Cedar Point, (MPN/lOO ml) 

total coliform concentration in the water of the 
Gulf of Mexico, (MPN/lOO ml) 

total coliform cell loading concentration at 
Mobile River mouth, (MPN/lOO ml) 

total coliform cell loading concentration at 
Tensaw River mouth, (MPN/lOO ml) 

depth of Hater in a cell, (ft) 

x- and y-direction grid size (=6561,68 ft) 

time increment (=240 sec) 

x-component dispersion coefficient, (ft^/sec) 
y-component dispersion coefficient, (ft^/sec) 
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Variable 

I 

IB 

IBNDL 

IBNDR 

IFLD 

IFRM 

IKK 

IPRNT 

I QUIT 
IRCB 

IREP 

ISTRT 

11 

12 

13 

14 
J 
K 

KIK 

NC 

NUM 


Name 


Definition 


counter 

number of rows possessing boundary cells 

cell numbers of the left hand side boundary cells 

cell numbers of the right hand side boundary cells 

specification used in variable format to designate 
position of decimal point 

specification used in variable format to designate 
start-printing position, number of variables, 
spacing, etc. 

counter 

specification used in variable format as printer 
control code 

number of the cell to stop printing in that row 

number of cells in each of the two ideal river 
channels (=10) 

number of cells in each row where results are 
to be printed 

number of the cell to start printing in that row 

counter 

counter 

counter 

counter 

counter 

counter 

counter 

number of cells in the grid system (=798) 
spacing to be indented before printing begins 
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Variable Name 

Description 

TCALP 

total coliform cell loading concentration at 
Alabama Port, (MPN/lOO ml) 

TGBSM 

total coliform cell loading concentration at 
the mouth of Bon Secour River, (MPN/lOO ml) 

TCDRM 

total coliform cell loading concentration at 
the mouth of Dog River, (MPN/lOO ml) 

TCHK 

number of tidal cycles elapsed after which 
computation is stopped 

TDEL 

time expended in computation (sec) 

TDLE 

number of tidal cycles elapsed (=25 hr. each 
period) 

TEMPF 

bay water temperature (°F) 

TEMPO 

bay water temperature (°C) 

TIM 

maximum time for run (=4800 hr. ) 

TPRNT 

time interval between printing of results 
(=360,000 sec) 

ULNT . 

x- component net velocity over one tidal cycle, 
(ft/sec) 

VLNT 

y- component net velocity over one tidal cycle, 
(ft/sec) 



A3 Program Listings 


(i) "CREATE NCSTM"(Storing of the Main Program of the NCSTM 

in the UNIVAC-1110 Computer System) 

(1) Card Arrangement (Control Cards) 

©RUN , A/TPC HUAAN , (account no . ) , LIU , 2 , J0/0 
@ASG,CP HUAAN ,F 
©FOR, IS HUAAN. NONC 

Card Deck of the Main Program 

©MAP, I , HUAAN. NONC 
IN HUAAN. NONC 
©PREP HUAAN . 

©PACK HUAAN. 

@PRT,I HUAAN. 

@@ 

Blank Card 

(2) Listing of the Main Program 
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NO w C ItalO 

C < :i /'j ) -COl.1 
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(il) "RUN NCSTM” (Gall the Main, Input with Data, and Compute.) 


(l) Card Arrangement (Control Cards) 

@RUN,A/TPC HUAAN, (account No.), LIU,2,30/0 
@ASG,A HUAAN 
@ASG,T 3 

©FOR , SU HUAAN . NONC , . NONC 


Updating Cards 

©PREP HUAAN. 

©PACK HUAAN. 

©MAP, SI , HUAAN. NONGXQT 
IN HUAAN. NONC 
©XQT HUAAN. BODXQT 

Card deck for the Input Data 

©ADD , P APRIL. PDATA 
©ADD , P HUAAN. IIDAT 
@ELT,IL HUAAN. IIDAT 
©ADD , P 3 
©PACK HUAAN 
@@ 


Blank Card 


* This control card serves to pick up the results obtained by 
the previous run as initial conditions. When the NCSTM is 
run for the first time, leave this card out. After the first 
run is completed, put this card in and insert the FORTRAN 
statement 

READ(5,8l) (C(I),I=1,NC) 

in "between statements HON Cl 050 and N0HC1 060 of the main 
program of the NCSTM as listed in the preceding pages. 


(2) Listing of Typical Input Data 




TVPICAL IhPUT DAT * FOR THF NCSTM FOR MOBILE BAY 
TYPICAL DATA CArlbS Or THE MCSTM FOR MOBILE BAY 

CAKO 1 IB*nC>IrCB 

36 . 7*6 10 


. Cnfi'J 


t>ah i • urt 


C«R;j 


uo. 


Ct.rVJ 


d OELS»u&L7.COLKfCOLT»TlM»Tt»MNT.COL®»COLC 
2^u . yOu4 » f>0 ft* 4u00 » 360000 • 

3 TCBSHt TCALPf TCDRM 

1 100 • 900 « 

v Tr.MPF 


7<» • 1 



LmkD S— liiHyL 


. 2 

23 

44 by 

(■9 109 

129 

149 

i7n 

1Q2 

213 

235 

256 

277 

29° 

316 

340 

361 

382 

403 

■' 7> 


407 489 

610 531 

yS2 

S74 

59 b 

61 6 

o37 

656 

679 

700 

721 

742 

763 

764 



4 

■ ».<0 

7-o 

1 AOijiK 















1 1 

yd 

S3 7u 

99 124 

145 

166 

160 

206 

226 

247 

266 

286 

306 

327 

348 

368 

389 

411 


4 64 

474 406 

016 637 

5SH 

679 

60 ft 

621 

642 

663 

604 

705 

726 

747 

7ofi 

789 



c 

kkD 

9-10 

Nj-j.v; 















6 

6 

«‘i n 

2b 2 0 

14 

rt 

6 

1 4 

14 

20 

20 

20 

20 

20 

20 

20 

60 

20 

?6 

26 

26 32 

32 32 

32 

3H 

3B 

36 

36 

35 

38 

38 

38 

38 

38 

38 



c 

. uRO 

11-12 

IrttP 















1 0 

10 

10 6 

11 la 

17 

ia 

IB 

15 

14 

13 

11 

10 

6 

9 

9 

8 

0 

9 

« 

9 

7 

7 7 

7 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 
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TypICmL iJvPUT DATA FOR THft NCbTM FOR MOBILE RAY 
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37 34 39 28 17 30 lo ?7 38 39 40 40 3? 30 32 33 37 39 35*37 

38 3^ 2j. 9 9 20 30 23 30 38 4 0 40 39 33 26 23 30 33 33 35 34 

31 29 o 11 9 23 30 

5 6 2o 20 10 ? I 1 3 b 5 

n 4 £> 9 14 l ? 14 12 3 7 7 9 7 ft 4 5 4 

. 4 <> 1» 12 15 12 12 12 12 12 10 9 7. ft a 6 5 

lu 10 11 12. 12 12 U 1U 14 10 11 11 10 7 9 0 B 5 

a rt 10 12 11 11 10 10 10 1G 11 HU 79 84 

7 10 11 10 Li 10 10 10 10 10 12 12 1? 8 3 

3 10 11 11 12 11 11 10 11 10 lft 9 7 3 

7 1 0 10 11 15 11 10 1 0 10 9 a‘5 1 

7 9 1 li Id ]0 il 11 l U 1 i) 8 2 

5 9 10 li 10 11 11 10 72 
u 10 lij v 12 1111 lit 5 

3 10 1 U 9 li 11 10 lit 4 

4 lo'iL) 9 1.1 12 9 7 4 

6 10 10 10 li 10 9 4 

8 9 9 Hi 11 10 n 6 

3 7 6 11 li 11 11 4 5 

7 8 9 19 10 U ID 9 

7 ft il 1C 10 lft 910 6 

6 7 ri lu li 11 10 7 

• O 1 0 i' 1 0 9 < i ft 

2 15 10 10 9 7 5 

2 710*ft:>10 

2 7 3 0 6 4 3 

u 2 4 ,3 1 6 

3 1 1 1 1 ft 

to 1 1 X 1 10 

lo 1 1 1 1 10 

15 1 1 1 1 10 

10 i L 1 l li) 

lu 1 1 1 1 10 

lo i i l i io 

to 1111 10 . 

10 1 1 • 1 1 10 

10 1 1 1 1 10 


*791 
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(ill) Computing Scheme 

The card sequence of the ’’CREATE NCSTM" was used to catalogue 
the main program of the NCSTM into the mass storage of the UNI VAC- 
1110 computer eye tern. The card sequence of the ’’RUN NCSTM" was used 
to call the main program, feed in data, and calculate the total coil- 
form distribution within Mobile Bay. Once the main program is cata- 
logued, only the "RUN NCSTM" deck is needed to effect computation. 
Changes In cell loading concentrations of total coliform and change 
of temperature for different months can be made by altering values of 
the data on data cards no. 2, 3» and 4 (see Appendix A4). Statements 
in the main program can be corrected, deleted, or the whole main 
program can be deleted from the mass storage by using appropriate 
control cards . 


(IV) Listing of Typical Output of the NCSTM 



R'J.: TOTkL COL .FORM FOR PARAMETRIC STUDY U 

TGTni. COt-iFORN iiACicRlA IN MPN PER 10 0 MI LI TER 
T- 12 ,'jOCYCElS uEi. f= 240. UO 








16000. 

3422. 

1522. 

492. 





6157. 

7 656 . 

3194. 

1251. 

456. 





4501 . 

5935. 

2496. 

1014. 

411. 





3369 . 

4544. 

2017. 

652. 

386. 





2590. 

3403. 

1664, 

729. 

354. 




1195. 

1947, 

2463. 

1365. 

630. 

339. 




450 , 

1440 , 

1769. 

1029. 

545 . 

313. 




440. 

1065. 

1193. 

771. 

456. 

274. 



374. 

434. 

723. 

832. 

602. 

370, 

235. 



365 . 

362 . 

522 . 

620. 

461. 

302. 

198, 



390. 

329 . 

4 07. 

469. 

329. 

249. 

lti5 • 



464 . 

295 . 

326. 

356, 

314. 

205. 

137. 



j 15 . 

273. 

207. 

27b, 

233. 

169. 

115. 



Vi 5 . 

261, 

223. 

221 . 

202. 

140. 

96. 



Iij 5 . 

251. 

I9i.. ieo. 

161. 

116. 

61 . 



?i-7. 

230. 

.165, 

149, 

127. 

96. 

69. 



'J i 3 , 

205, 

192. 

123. 

1C1. 

79. 

59, 


Oil . 

5-5. 

ift 3. 

123. 

93. 

6K 

65. 

51. 


c 1 5 , 


1 5^ « 

104 . 

71 . 

63. 

54. 

44. 

107. 

132. 


121. 

6b. 

57, 

50. 

44. 

37. 

0. 

50 . 

■35 . 

69. 

70. 

46, 

41. 

37. 

31. 


0 5 . 

69 . 

65 . 

56 , 

37. 

35. 

32. 

26. 



So . 

5 1 

44. 

30. 

30. 

26. 

25. 




3o . 

34 , 

25. 

25. 

25. 

24. 

O • 

5 » 

10. 

21. 

21 . 

19. 




♦ 

4 . 

5 , 

a. 

11, 

12. 

4. 

2. 

1 . 

Cj « 

i . 

2 « 

4 . 

6. 

7, 

3. 

1, 

0, 

0 , 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


204 . 200 . 

222. 1B6. 

220 . 173 . 

21b. 15S. 

212 . 152 . 

£ 09 . 150 , 

200. 152. 136. 

164. 147. 

162. 136. 

140. 

116. 

99. 34. 

S3 . 6b • 

69. 53. 

S3. 43. 37. 


50. 

37. 

30. 

25. 






44. 

33. 

25. 

20. 

15. 

7 . 




39. 

29. 

22. 

17. 

' 12. 

a. 




35. 

27. 

20. 

15. 

11. 

a. 

6 * 



30. 

24. 

19. 

14. 

10. 

a. 

6 • 

6 . 


26, 

22. 

17, 

13. 

10.' 

a. 

H. 

10. 

15 

23. 

19. 

16. 

12. 

10. 

a. 

10. 

17. 

40 

21. 

17. 

12. 

10. 

a. 

7 . 

b« 

6. 

11 

21. 

17. 

a. 

6. 

5. 






0 . 

6 . 


0 


H 

On 

ON 
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A4 input data specifications 




Card 

No. 

Column 

No. 

Variable* 

Name 

Input 

Unit 

Input 

Format 

i 

1-5 

IB 

— 

Integer 


6-10 

NC 

— 

Integer 


11-15 

IRCB 

— 

Integer 

2 

1-10 

DELS 

ft. 

Real 


11-20 

DELT 

sec. 

Real 


21-30 

C0LM 

MPN/lOOml 

Real 


31-40 

COLT 

MPN/lOOml 

Real 


41-50 

TIM 

hr. 

Real 


51-60 

TPRNT 

sec. 

Real 


6l-?0 

C0LG 

MPN/lOOml 

Real 


71-80 

C0LC 

MPN/lOOml 

Real 

3 

1-10 

TCBSM 

MPN/lOOml 

Real 


11-20 

TCALP 

MPN/lOOml 

Real 


21-30 

TCDRM 

MPN/lOOml 

Real 

4 

1-21 

TEMPF 

°F 

Real 

5-6 

every 4 columns 

IBNDL 

— 

Integer 

7-8 

every 4 columns 

IBNDR 

— 

Integer 

9-10 

every 4 columns 

NUM 

— 

Integer 

11-12 

every 4 columns 

IREP 

— 

Integer 

13-50 

every 3 columns 

D(I) 

ft. 

Real 


* For description of the variables see Appendix A2. 



APPENDIX B 


EXAMPLES OF PAW DATA USED IN 
THE NCSTM 


Bl Total Coliform Bacteria 

B2 River Flow Rates 

B3 Wind Conditions and Temperatures 
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B1 EXAMPLES OF 
RAW DATA OF 

TOTAL COLIFOEM CONCENTRATION 
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_The Significance of JSC Positive Organisms _in Gulf Shellfish Growing 
Waters - H. S. Hosty, Alabama Health Department, Montgomery, Alabama. 

The examination of Mobile Bay was divided into three separate 
areas during the whole course of this study. Phase 1 . was to deter- 
mine the significance of conforms versus fecal coliforms as an indi- 
cator of pollution. Phase 2 . was an attempt to use pathogenic E.coli 
as an indicator of human pollution and Phase 3 . was a comparison of the 
sanitary quality of oysters harvested and shucked in the laboratory as 
compared to those harvested and shucked in individual plants. 

In Phase 1 all of the procedures were those recommended by the 
Bacteriological Examination of Seaweed and Shellfish, third edition 

/ 

1962. This investigation involved the weekly testing of A3 stations for 
two years or approximately 4500 coliform examinations by the three- 
tube test. Additionally, all positive lactose tubes were reinoculated 
into EC media and incubated in a water bath at 44.5° = 0.2° for 24-48 
hours. For sometime, after incubation, all EC tubes positive or 
negative were streaked to eosin methylene blue plates. After incubation, 
one colony conforming to the accepted morphology of E.coli, or if no 
such were present, a colony which came closest to being typical, was 
picked and inoculated into a lactose tube. From this tube the EC test 
was repeated along with IMVIC determination. It soon became apparent 
that plating of EC negative tubes was not fruitful so the routine was 
altered, plating only positive EC specimens. In all some 20,000 IMVIC 
and repeat EC determinations were performed. No tests were performed 
on oysters but this is now under study. 
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The Mobile Bay area is roughly some 30 to 35 miles long and approxi- 
mately 25 miles at its widest point. The upper part of the bay shows 
extreme pollution which has no effect on the Cedar Point growing areas 
under normal conditions. During the rainy season this immunity to con- 
tamination may abruptly change. It may be added here that last February 
and March Alabama experienced the worst floods in recorded history. 

Figure 1 shows a general picture of the station locations as well 
as the normal current flow and flood pattern. A few of the upper 
stations, including the Alabama and Tombigbee stations, are not shown 
in this chart. By studying the various stations it was apparent that 
some areas could be grouped rather than considered independently. 
Alphabetically, therefore, grouped or single stations, were, in their 
descending order down the bay, starting with the Alabama and Tombigbee 
Rivers, designated as Group A, B, etc. In general we shall only discuss 
Group H, J, L and M. The station or stations represented by this group- 
ing may be identified in the table showing probability percentiles. 

Figure 2, Station 31, adequately shows the high pollution usually 
present. The overall trends suggest diluting out of EC gas positive 
organisms derived from the fresh sewage discharges by fresh water during 
periods of increased river flow. This is apparent beginning in late 
October, with progressively increasing separation of the coliform and 
EC lines throughout the flood stage which persisted until late May. 

When normal river flows reoccurred „ there was a return to a similar 
pattern of the coliform and EC MPN's. 

Station 48 (Figure 3) is included because it represents an area 
affected by run-off which reached a peak towards the end of March. The 



1?2 


rythmical pattern of coliforra and EC MPN separations were disrupted 
when salinity decreased. The marked increase in run-off resulted in a 
doubled coliform count but the EC MPN was only slightly changed. This 
continues until recovery of normal salinities is apparent. 

Some 35 miles down the bay from Station 31, Station 88 (Figure 4) 
has been used as an index area. During normal salinity content in the 
bay the coliform count, though fluctuating, reaches sizeable peaks while 
the EC count remains less than three. Around the 20th of March, with 
salinity, pollution occurred and there is an immediate rise in 
EC MPN s» The pattern of recovery shows gradual return to the usual 
pattern of coliform-EC relationships. 

Station 119 (Figure 5) is directly over the Cedar Point oyster 
area and illustrates the response of the EC test to pollution. The con- 
current drop in salinity and rise in EC MpN f s is dramatically evident. 

The histogram (Figure 6) summarizes in percentage the confirmed 
EC positives as compared to coliform positive tubes and shows clearly 
the response of EC media to pollution as compared to the accepted 
coliform test. During period 1 at Group G (Station 20) immediately 
below the discharge of fresh sewage, the EC test was 83 percent positive 
compared to the coliform test and there was a gradual decline until 
Group M (oyster bed area) was reached where only 9.8 percent of the 
coliform positive tubes were EC positive. During the period 2, flood 
stage, no change in percent EC positiveness was apparent in area G. 

Some dilution by river water caused the reduction in percent EC positives 
at groups H and J. The most striking difference occurred in the areas 
immediately above and over the oyster beds. Sharp increases in percent 
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EC gas positives occurred. These can be explained by the influx of 
fresh water that received pollution in the Mobile area and, as shown in 
Figure 5, reached the oyster beds during the flood stage. The fact that 
travel time was considerably less during the flood stage with the 
associated reduction in dieoff rate might also be a factor which 
influenced the increases in percent EC gas positiveness. It is also 
significant that when river discharges returned to normal low rates 
(recovery period), the progressive decrease in percent EC positiveness, 
associated with distance from pollution, was demonstrated. 

Table 1 shows these same grouped stations as coliform and EC 
percentile probabilities. Again, as one descends the bay, improvement 
is noted over the oyster beds (Group M) the 50 and 90 percentiles fall 
well within accepted norms during the normal and recovery periods but 
completely outside the range during the flood stage. 

From 8,400 positive E.C. tests incubated for 24 to 48 hours the 
following was recovered: 24 hour period 88.4 percent types 1 and 11 

E. coli strains, 7.9 percent irregular Vi's and 1.3 percent as other. 

By contrast incubation for an additional 24 hours resulted in the 
recovery of only 0.2 percent more E. coli but an additional two percent 
were irregular Vi's and others. Should this trend continue it seems to 
be that we are lowering the specificity of the test without increasing 
the sensitivity by the additional 24 hour incubation. 

Some 600 strains of types 1 and 11 E. coli and 181 isolates of A. 
aerogenes were subjected to the E.C. test run at 44.5 and 45 C. Raising 
the temperature one half degree eliminated 60 percent of the aerogenes 
isolates while less than one percent of the E. coli strains were lost. 
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Figure 1 . Coliform Sampling Station Locations in Mobile Bay. 
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RELATION of COLIFORM POSITIVE TUBES to CONFIRMED 
EC and IMViC TYPES by PERIODS 
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B2 EXAMPLES OP 


RAW DATA OP 


RIVER PLOW RATE 



ORIGINAL PAGE IS 
OP POOR QUAUTSi 


Table B2.1 

2-4705 MOBILE HI VO- NEAR MOUNT VERNON ALA 
LOCATION— LAT 31 06 50 LONG K7 5ft 05 IN SE 1/4 SEC 41 T2N RlE AT BOAT PIER OH 
hCST (>«, OF LAKE DAVID .5 mI US FR LAKE OUTLET°TO MOBILE RV 2.5mI NE OF HT VERNON 

+ AT MtLt 41.3 FR MOBILE /DRAINAGE AREA 43000 /HECOROS AVAILABLE OCT 1953 

.TO SEPT 1054 / GAGE STAGE RECOVER DATUM AHOUT 2FT BELOW MSL UY COMPARATIVE 

Gil READINGS AUX GAGE AT ALA ST DOCKS 

Thru ivg7 mount vernon flow =i claiuorne flow + hy 43 flow! * 1.05 

AFTER 196 7 MOUNT VERNON FLOW =t CLAlfiORNE FLOW ♦ COFFEEVILLE FLOW) *1.07 
monthly and yearly discharge in cubic Piet per second 
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Table 32*2 

2-9705 MOBILE IUVd NEAR MOUNT VERNON ALA 
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TO 5£PT 1954 / GAGE— -STAGE RLCOiu-LR DATUM ABOUT 2FT BELOw M5L bY COMPARATIVE 
GH REAOINGS AUX GAGE AT ALA ST DOCKS 

THRU 1907 MOUNT VERNON FLO W st CLAlyORWE FLOW + HY 43 FLOW) * 1.05 
AFTER I9ia7 MOUNT VERNON FLOW =1 CLAljiORNE FLOW + COFFEEVILLE FLOW) * 1.07 
A MONTHLY AND YEARLY DISCHARGE IN CUBIC FLET PER SECOND 

OCTOBER NOVEMBER DECEMBER JANUARY FEBRUARY MARCH APRIL 
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CLIMATOLOGICAL DATA 
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Fig. B 3 .I Bimonthly surface isothermal maps of Mobile Bay and Mississippi Sound 
Alabama, Combined data from 1963-64, 1965-66 (McPhearson, 1970) and January 1968 
through March, 1969. 
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Fig* B3.2 Bimonthly bottom isothermal maps of Mobile Bay and Mississippi Sound, 
Alabama. Combined data from 19G3-64, 1965-60 (McPhearson, 1970) and January, 1968 
through March, 1969. 



